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1. INTRODUCTION 


Over the life span of this NASA-VAS cooperative program, work has been 
divided between analysis/forecast model development and evaluation of the 
impact of satellite data in mesoscale numerical weather prediction (NWP) and 
also development of the Multispectral Atmospheric Mapping Sensor (MAMS) and 
related research. 

The modelling effort has seen the CIMSS Synoptic Scale Model (SSM) 
progress from a relatively basic analysis/forecast system at the inception of 
the program to a package which now includes such features as nonlinear 
vertical mode initialization, comprehensive Planetary Boundary Layer (PBL) 
physics and the core of a fully four -dimensional data assimilation package 
which will be expanded on for subsequent NASA sponsored research. 

Investigations on the use of satellite information have been rewarding in an 
educational sense, but sometimes frustrating in practical application, in that 
the data have not had the impact on mesoscale prediction which the research 
community had hoped for. In retrospect, the analysis and forecast problems 
associated with characteristics of the data could have been understood earlier 
by an amalgamation of expertise from the data producers and analysis - 
forecasting groups. It has been this NASA program and similar efforts which 
have inspired the community to consolidate the experience of these groups into 
a common ground of expertise which will enable more timely and beneficial use 
of the information from future remote sensing data sources. 

The MAMS effort has produced a calibrated visible and infrared sensor 
that produces imagery at high spatial resolution. The MAMS has been developed 
in order to study small scale atmospheric moisture variability, to monitor and 
classify clouds, and to investigate the role of surface characteristics in the 
production of clouds, precipitation, and severe storms. The NASA-VAS 
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cooperative program has been the starting point for the testing of this 
aircraft instrument; the design of a future space borne sensor in low earth or 
geostationary orbit with similar monitoring capabilities in planned with EOS 
and beyond. The effort funded here has demonstrated good quality data from 
the MAMS from intercomparison with other instruments. 

2. SATELLITE DATA ASSIMILATION RESEARCH 

With the satellite era has come the first mesoscale data base available 
to the meteorological community on a routine basis. While these satellite 
systems have made sizable contributions to the understanding of global -scale 
processes, contributions to mesoscale meteorology have been slower to evolve. 
Both case study work and observing system simulation experiments (OSSES) echo 
the same theme. On the mesoscale, traditional methods of looking at his data 
base must be modified and not only the data errors and sampling 
characteristics, but also the physical meaning of the measurements must be 
well -understood. Though research approaches vary between investigators, there 
is a common underlying theme to satellite data assimilation. It is evident 
that, almost without exception, the techniques have been designed to access 
the information content of the data while circumventing one or more of the 
problematic aspects of the observing system. 

Table 1 contains a brief summary of qualitative characteristics of the 
TOVS, VAS and cloud and water vapor drift wind (CDWVW) measurements which have 
somewhat subjectively divided into "advantages" and "disadvantages". Most of 
these are well-known, but will provide a convenient context in which to 
discuss some of the methods used to analyze and assimilate the data. 
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Table 1. 


Characteristics of Satellite Data 


TOVS 

VAS 

CDWVW 

Advantages • 

moderate/high 

horizontal 

resolution 

high horizontal 
resolution 

high horizontal 
resolution 
where tracers 
exist 

regular coverage 
at set times 
depending on 
orbits 

user-selected 
frequent coverage 
in clear areas 

frequent coverage 
if rapid-scan 
mode used 

spatially 

correlated 

errors-good 

gradient 

information 

spatially 

correlated 

errors-good 

gradient 

information 


Disadvantages : 
low vertical 
resolution 

low vertical 
resolution 

level assignment 
uncertainties 

different 

error 

characteristics 
for clear, 
partly and 
cloudy retrievals 

gaps in coverage 
in cloudy areas 

single level 
only 

spatially 

correlated 

errors 

(biases) 

spatially 
correlated 
errors 
(biases ) 



A. TEMPERATURE AND MOISTURE RETRIEVALS 


On the global scale, TOVS and VAS temperature and moisture retrievals may 
be treated similar to their radiosonde counterparts and show generally 
positive impact on global forecasts , “especially in areas of low conventional 
data density such as the southern hemisphere. In these data-poor regions, the 
satellite data represent an improvement in the designation of atmospheric 
state over that achieved by conventional synoptic information and lower 
resolution forecast models alone. However, most of the operational global 
scale assimilation schemes use statistical interpolation methodologies in 
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which satellite data are allocated greater observational error and differing 
covariance functions than radiosonde data and also utilize data quality 
algorithms so that a high quality observation cannot be overwhelmed or 
rejected by a larger number of satellite observations of lesser quality* 

On the mesoscale, however, and especially in regions of relatively high 
density synoptic data and high resolution forecasts such as the continental 
United States or Europe, the characteristics of satellite soundings can lead 
to a degradation of forecast quality when they are assimilated in a standard 
scalar fashion. Our research efforts on the mesoscale have concentrated on 
methods to access the information content of the soundings or radiances while 
circumventing their problematic features. 

Probably the most detrimental feature of VAS and TOVS soundings to 
mesoscale assimilation work is their biases, which can be large and situation 
dependent. The biases arise from a number of sources. The limited vertical 
resolution of these satellite instruments (optimistically five levels of 
vertical information) immediately dictates that atmospheric features with 
small vertical scales may be missed or only accounted for in a vertically 
integrated manner. A second source of bias error comes from the retrieval 
algorithms themselves. Implicit or explicit in these algorithms, whether they 
be statistical or physical, is a comparison of measured radiances to "forward" 
radiances calculated from the background atmosphere for the retrieval. Errors 
in the knowledge of atmospheric transmission, instrument calibration error, 
errors in the assignment of surface emissivity, undetected sub- field-of -view 
clouds for infrared instruments and other errors may all enter and degrade the 
retrieval process. Even when the background field for retrieval is 
essentially "perfect" (a high resolution analysis from raob data) , we have 
shown that errors associated with the retrieval process may significantly 
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alter the satellite soundings to the point where they cannot be used for 
mesoscale NWP. 

Our research to mitigate this problematic aspect of satellite soundings 
has two facets. At the assimilation end, we have concentrated on developing 
methods for the assimilation of the horizontal gradients of temperature and 
moisture from VAS atmospheric soundings using variational techniques . 

Diak (1987) investigated the assimilation of VAS information by both 
scalar and gradient methods for the AVE/VAS case study day of 6 March 1982. 

On this day, there was a mesoscale network of raob observations in the Texas 
region for verification of analyses and forecasts. It was a day of large 
temporal changes in the region from a fast-moving cold frontal system and also 
a day in which large VAS data biased existed. Figure 1 shows 1500 UTC 
analyses of 500rab height errors for VAS data versus special network raobs 
using the VAS information as scalar (1A) and gradient (IB) information, 
respectively. As shown, the gradient approach significantly reduced error 
over the region compared to the scalar analysis. 

A forecast made from the gradient (1500-0000 UTC) was superior to the 
corresponding scalar forecast and similar in quality of primary variables to a 
forecast (1200-0000 UTC) initialized with only synoptic data. Improvements 
were noted in vertical motion and precipitation patterns in the gradient 
versus the synoptic forecast. More details on this experiment can be found in 
the enclosed reprint. 

While gradient methods mitigate some of the problems of biased data, it 
can be seen in principle that they are not a panacea. Gradient methodologies 
essentially solve a boundary value problem; that is, while gradients within a 
data region may be well-depicted, the absolute values of an observable in the 
region will depend on the invariant values of the background field on the data 
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boundaries. Our research has indicated that the gradient methodology performs 
best when satellite gradients are combined with measurements from raobs to 
"lock in" the absolute values of an observable at the synoptic locations 

Figure 2 shows an example of where this technique may be useful in the 
western Atlantic Ocean region. Figure 2a illustrates the problem of trying to 
produce a realistic analysis of moisture with the limited conventional data 
available in the western Atlantic. As shown, the result of limited data in 
our analysis system and others such as the LFM and RAFTS is often the "Bermuda 
Bulls Eye" where the Bermuda raob and the background field for analysis (in 
this case, the NMC Global forecast) are at odds. Figure 2b is an analysis 
made from a combination of satellite gradient and raob absolute value 
information. The satellite gradients have filled in the data voids between 
synoptic locations and improved the coherence of the moisture analysis in this 
case . 

In addition to gradient analysis methods, some research progress had been 
made in mitigating data biases through retrieval techniques which rely on 
measurements of the time change of radiances. The method acknowledges that in 
areas of dense synoptic station coverage the best depiction of the atmosphere 
is at synoptic times. The retrieval problem can then be formulated as the 
time change of atmosphere (radiances) from this well-known initial state, that 
is : 


“ ®F + 

Bp - radiance for atmospheric retrieval 

Bp - "forward” calculation of radiance at synoptic time 
AB^ - measured change in radiance (data time-synoptic time) 

Very important in this retrieval methodology is that bias errors from 
instrument sources, atmospheric transmissivity errors, surface emissivity 
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errors and other sources, drop out of the retrieval process in a first order 
fashion. 

At this writing, a retrieval algorithm using the time change of radiances 
has been written and required only minor modifications to the standard CIMSS 
physical retrieval method to develop. It is planned that investigation of 
this technique will proceed under other existing NASA/CIMSS programs. 

Some positive preliminary results have been drawn from a conceptually 
very similar experiment (synthetic radiances and soundings) using the time 
change of the VAS retrievals themselves. Figure 3a shows an initial 500mb 
temperature state analyzed from raob special network observations at 1500 UTC 
6 March 1982. Figures 3b and 3c are two satellite analyses of VAS retrievals 
of 500mb temperature at 1800 UTC on the same day. The first is an analysis of 
VAS retrievals made in the standard manner, while the second analysis was made 
using the 1500 UTC initial state (2a) augmented by the time change of the VAS 
retrievals 1800 minus 1500 UTC. All retrievals were made using an LFM first 
guess. Comparing these two analyses to the 1800 UTC verification (3d), 
another special network raob analysis, it can be seen that the temporal scheme 
was superior in maintaining small-scale features and the strengths of the 
temperature gradients in the region. The time-changes (4b) show reasonable 
correspondence to those from the special network analyses (3a) and importantly 
add information to the time changes depicted in the LFM guess fields which 
were used to make these soundings (4c) . 

It can be seen that the time change of atmospheric radiances can be a 
useful tool, not only in dealing with satellite data biases, but in the more 
general framework of understanding the coupling of thermodynamic and radiative 
processes in the atmosphere. 
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Perhaps the poorest understood element of satellite soundings in an 
assimilation context is the effects of their spatial and temporal sampling 
characteristics . Polar orbiting instruments sample atmosphere radiances at a 
given location only once or twice per day, while geostationary infrared 
instruments are subject to large data voids caused by clouds. On the global 
scale it is standard procedure to bin their observations into a designated 
time window of about six hours for assimilation. One mesoscale experiment at 
the CIMSS has shown that as little as a two hour difference in time between 
data from raob and TOVS sources noticeably degraded the quality of an analysis 
of a fast moving cyclone made from the combination of the two data types. 

We have achieved some success using complementary data sources to fill in 
the data voids of infrared sounder. LeMarshall et al . (1985) and Diak et al. 

(1986) both describe the complementary nature of the spatial distribution of 
VAS soundings and cloud and water vapor drift winds. In these studies, the 
CDWVW were used as indicators of geopotential gradient and helped describe the 
mass field in sounding data voids. Unfortunately (as will be described in the 
next section), the single-level nature of these winds limits their routinely 
demonstrating positive impact on numerical analyses and forecasts. 

B. CLOUD DRIFT AND WATER VAPOR WIND DATA 

There are two problems which can be immediately identified with the use 
of CDWVW data in analyses and forecasting. The first is one of height 
assignment. Given that the cloud or water vapor tracer has some depth, to 
what level should the wind vector associated with the tracer be assigned. 
Multispectral techniques, more sophisticated that simply assigning low middle 
and high levels, are in common use, but for the time being at least many 
assimilation schemes make some allowance for the uncertainty in level 
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assignment by allocating a greater observational error to the data in 
statistical interpolation and by using appropriate data rejection algorithms, 
particularly if the observation is near a jet stream. 

The second problem is that these data are at single levels and cannot by 
themselves provide a three-dimensional wind analysis, foreshadowing the 
possibility that they will not be well accepted by a forecast model. The 
variational analysis scheme which we employ does spread the influence of these 
winds in the vertical, however, not as effectively as an explicitly three- 
dimensional optimal interpolation procedure. 

A feature of the CDWVW data, as pointed out by LeMarshall et al. (1985) 
and Diak et al . (1986), and mentioned in the previous section is that they 
frequently complement the coverage of infrared satellite soundings Figures 
5a and 5b show, respectively, the VAS and cloud and water vapor drift wind 
distribution at 1600 UTC on the AVE/VAS case study day of 26 April 1982. On 
this day, a small closed low developed at 500 mb in the southern Illinois 
region between synoptic stations and synoptic times (evidenced in the cloud 
circulations). The feature did not develop in a 12-hour mesoscale model run 
from a synoptic analysis made at 1200 UTC 26 April. Figures 5c and 5d are two 
500mb analyses made at 1600 UTC using this forecast as a guess. The first 
(5c) is an analysis which used only the VAS soundings (as gradient 
information, see prior section) as input, while the second (5d) contains both 
VAS and WVCDW data. As shown, the CDWVW information was instrumental in the 
analysis of this feature. Three similar pairs of analyses for later times 
(with and without CDWVW data) yielded the same character of results. While 
the quality of results for this case is somewhat fortuitous owing to the good 
horizontal and vertical distribution of satellite wind data on April 26, it 
suggests that these data may be of at least periodic benefit when assimilated 
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in a consistent manner with other observations. LeMarshall et al . (1985) 

similarly found improvements in the forecast position of a tropical cyclone 
when both VAS and CDWVW data bases were used in the analysis. 

We expect that the forthcoming upgrade of our analysis system to a 
combination of optimal interpolation-variational blending will increase the 
frequency with which CDWVW data display positive analysis and forecast 
results . 

3. MULTISPECTRAL ATMOSPHERIC MAPPING SENSOR (MAMS) 

This final report will describe the characteristics of the instrument 
design that evolved from this study, present the results of some of the VAS 
and MAMS intercomparisons, and summarize the science data that has been 
gathered. This work is continuing in a follow on contract with NASA, so that 
a more complete final report will be forthcoming at the end of that contracted 
effort. 

A. MAMS INSTRUMENT CHARACTERISTICS 

The MAMS is a modification of NASA's Airborne Thematic Mapping (Daedalus, 
1982) that was suggested by UW. The MAMS has two major components, the scan 
head and the spectrometer. The scan head consists of the primary collecting 
telescope, a rotating scan mirror, a motor encoder assembly, and two 
controlled thermal sources. The field of view is defined by the aperture 
which is available in two different sizes: 2.5 mrad and 5.0 mrad. The 

spectrometer consists of optical elements which spectrally separate the 
polychromatic input energy, lenses which focus the separated energy onto 
sensors, detectors which convert optical energy to an electrical signal, and 
pre -amplifiers which condition the signal. The energy collected by the 
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spectrometer is separated into four distinct optical paths using dichroic 
filters; three paths for the infrared data (one for each spectral band) and 
one for the visible and near infrared data. The spectrometer detectors are 
silicon for the visible radiation and mercury cadmium telluride for the 
thermal infrared radiation. The visible detector is an eight element array 
which converts optical energy to electrical energy. The spectral response of 
each band is determined by the dispersion characteristics of the prism placed 
in the optical path and the location of the array in the dispersed beam. The 
thermal infrared detectors are each housed in a vacuum dewar which contains a 
cooled longwave filter which defines the particular spectral band. Each 
infrared detector has an associated pre-amplifier to process the electrical 
signal. The spectral bands are summarized in Table 2. The spectral response 
functions for the three infrared channels of the MAMS were specified to be 
similar to the 6.7, 11.2, and 12.7 micron channels on VAS . 

Table 2 

Range of MAMS Spectral Bands 

8 visible 3 infrared 

(microns) (microns) 


.42- .45 
.45- .52 
.52- .60 
.60- .62 
.63- .69 
.69- .75 
.75- .90 
.90-1.05 


6.20-6.85 (upper level water vapor) 
10.32-12.02 (window) 

12.20-12.56 (lower level water vapor) 


The MAMS produces high resolution imagery in these visible and infrared 
spectral bands from a high altitude aircraft at 20 km. The horizontal ground 
resolution of each individual field of view is 100 meters (50 if che smaller 
2.5 mrad optics are used). The total field of view for one scan is 86° which 
covers the width of the entire cross path of roughly 40 km. The scan rate is 
6.25 revolutions per second (12.50 when the 2.5 mrad optics are used) and 
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produces 67% overlap, 1-P, of each scan line with the previous one for nominal 
aircraft speed (740 km/hr). On each scan line, 716 samples are taken so that 
there is 52% overlap of neighboring pixels along a scan line (4% when the 2.5 
mrad optics are used) . These overlaps can be used to reduce the noise in the 
scene by spatial averaging. 

Radiometric calibration of the infrared data is accomplished by viewing 
two blackbody sources and by assuming the detector response is linear with 
respect to radiance (a fractional nonlinearity of less than 5x10'^ is 
observed) . Calibration occurs every scan line by viewing a warm and than a 
cold blackbody of known temperature. For a given spectral band, the radiance 
is determined from the temperature through the convolution of the normalized 
spectral response and the Planck function. The calibrated radiances are 
transformed back into temperatures using the inverse Planck relation with the 
nominal spectral band center wavelength. 

Prior to flight, the gain and offset for each spectral band are adjusted 
electronically in order to minimize the likelihood of saturation over hot 
scenes and to maximize the effective operating range. The 6.5 micron band is 
adjusted to cover 210 to 280°K, while the 11.2 and 12.3 micron bands operate 
from 220 to 340°K. An eight bit digitizer is used so that for the 6.5 micron 
band, the gain is approximately .3°K per count, and for the remaining surface 
viewing bands it is approximately .5°K per count Truncation errors are 
assumed to be roughly half these values. 

The infrared detectors on MAMS exhibit some low frequency noise so that 
samples taken at 52 microsecond intervals for the spin rate of 6.25 
revolutions per second (26 microsecond sampling occurs at 12.5 rps) experience 
noise contributions with significant correlation. An analysis of covariance 
of the noise with the 5 mrad aperture reveals that after averaging ten 
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consecutive samples, noise is reduced by a factor of .54, .52, and .37, 

respectively, in the three infrared spectral bands in order of increasing 
wavelength (uncorrelated noise would have a factor of 1/710 or .32). 

Inflight single sample noise values were determined from multiple samples 
over a uniform portion of the ocean. For a sample size of 100, the single 
sample noise values with the 5 mrad aperture were found to be .2, .1, and 

.4°K, respectively, for the spectral bands in order of increasing wavelength 
(for the 2.5 mrad aperture the values were .5, .2, and .9°K). 


B. DATA INTERCOMPARISON 

Radiances from the MAMS and VAS spectral responses were simulated for 
each band and sensor in a radiative transfer calculation. A representative 
atmospheric transmittance was computed using a band model and radiosonde 
temperature and moisture profiles. Table 3 shows that the MAMS detects 
radiances within one to two mW/ster/m /cm of the VAS detected radiances. 
The MAMS spectral responses are shifted to shorter wavelengths, hence less 
radiance is detected. This converts into brightness temperature differences 
of less than a degree Kelvin. 

Table 3 

Intercomparison of Simulated Data 



6 micron 

11 micron 

12 micron 

R(MAMS) 

7.86 

96.5 

109.5 

T(MAMS) 

254.8 

287.3 

286.0 

R(VAS) 

8.73 

97.5 

111.6 

T (VAS ) 

253.9 

287.3 

285.4 

at mams-vas 

0 . 9 °K 

0.0 

0.6 

radiances 

2 - 1 

in mW/ster/m /cm , 

temperatures in 

°K 
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Radiances from the MAMS and VAS were collected simultaneously over a 
relatively isothermal area of ocean off the California coast on 18 May 1985. 
Table 4 presents the comparison. After compensating for the inherent spectral 
differences of the two instruments, the MAMS data is within .5, 2.5, and 1.5°K 
of the VAS data, respectively, for the spectral bands in order of increasing 
wavelength. The observed split window channel differences (11 micron less 12 
micron) are comparable, 2.5°K for MAMS and 2.1°K for VAS and thus showing 
similar low level moisture concentrations. The observed water vapor channel 
data (6 micron) show MAMS detecting upper level moisture concentrations very 
much like VAS. with MAMS showing somewhat less attenuation than VAS. The 
comparisons are quite reasonable, since the spatial resolution of each sensor 
is very different (VAS at 7 km and MAMS at .05 km at nadir), the navigation of 
the MAMS data is from visual siting of landmarks only, and non-nadir viewing 
Of the MAMS has not yet been taken into account. 

Table 4 

Intercomparison of Observed Data 


6 micron 

R(MAMS) 6.01 

T(MAMS) 247.1 

R(VAS) 6.80 

T(VAS) 246.5 

AT MAMS-VAS 0.6°K 


11 micron 

99.9 

289.5 

97.0 

287.1 

2.4 


12 micron 

111.1 

287.0 

111.1 
285.0 

2.0 


radiances in mW/ster/m^/cm ^ , temperatures in °K 


Imagery from the MAMS and VAS were compared from flights over Oklahoma 
and Kansas on 12 May 1985. The goal of the intercomparison is to verify small 
scale features observed in the VAS water vapor images with the high resolution 
MAMS data. These water vapor inhomogeneities in the VAS data have been used 
to obtain motion vectors when several images are considered in sequence. 
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While the imagery is often too amorphous to track winds reliably with a 
correlation technique, it has produced good wind fields by single pixel 
tracking. The MAMS data tends to verify that small scale gradients (down to 
several single pixels) in the VAS data are real. 

C. MAMS DATA ARCHIVE 

The MAMS has been flown on a NASA U2/ER2 high altitude aircraft for both 
engineering checkout flights and scientific missions. Flight dates and 
general information are listed in Table 5. While the January 1985 flights 
served primarily as engineering checkout of instrument and/or hardware 
problems, many interesting atmospheric and land features were observed with 
the instrument. In particular, mountain wave features in the lee of the 
Sierra Nevadas were observed in the water vapor imagery from 22 January (see 
Jedlovec paper on this subject). 

The MAMS flights during May and June 1985 were conducted to sample a 
variety of weather conditions and also corresponded to a period of special 
ground-based observations of the Pre- STORM field experiment. 

Several days are of particular interest. The first part of the 18 May 
flight contained ocean data, which provided an ideal situation to study the 
noise characteristics of the MAMS over a relatively stable background scene. 
Corresponding data from the VAS and AVHRR were available over this region at 
approximately the same time. Also, 20 and 24 May and 26 June are interesting 
because they sampled a humid summertime environment preceding thunderstorm 
activity. Flights on the 22 and 23 June were constructed such that four 
passes over the same area region occurred at one hour intervals, thereby 
providing high horizontal resolution with increased temporal frequency to 
monitor environmental changes . 
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The MAMS flights in January 1986 tested several new configurations for 
improving the signal - to-noise ratios. The speed of the scan head was halved 
(as well as the digitizing rate) and the 5.0 mrad aperture was thoroughly 
tested. From the data of these four flights, it was apparent that the signal - 
to-noise was substantially improved. 
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Fig. 1. Errors in 500 n 
and gradient assimilatio 
raob analysis 1500 UTC 6 



Is (gpm) for scalar assimilation (la) 
soundings versus special network 





FIG 2 

Fig. 2. First guess (radiosonde analysis) of 700 mb dewpoints and 
combination analysis resulting from variational blending of horizontal 
gradients from VAS data into first guess. 
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FIG 4 


Fig. 4. Time change of 500 mb temperature (K) 1800 minus 1500 UTC 
6 March 1982 from special network raobs (3a), VAS (3b), and LFM 
forecast (3c) . 


JR133MAL PAGE i'S 

OF POOR QUALny 



FIG 5 


Fie. 5(a). VAS sounding locations 1600 UTC 26 April 1982. (b) CDWVW 

locations 1600 UTC 26 April 1982. (c) Analysis of 500 mb geopotential 

(gpm) using VAS soundings 1600 UTC 26 April 1982. (d) Same as (c) , 

only using VAS plus CDWVW. 
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abstract 

Comparisons are made between analyses and forecasts which incorporate VAS geopotential data as either 
scalar or horizontal gradient information for a case study on the AVE/VAS day of 6 March 1982. On this day, 
incorporating the VAS information in analysis as a variational constraint on horizontal geopotential gradients 
significant mitigated the effects of large data biases which made VAS assimilation by standard scalar methods 
very difficult. A subsequent forecast made from the gradient assimilation was superior to one made from the 
standard analysis and of comparable quality in geopotentiais to a control forecast from synoptic data. Most 
impact was noted in the forecasts of vertical motion and precipitation in the gradient versus this control simulation. 


1. Introduction 

During the past several years, the Cooperative In- 
stitute of Meteorological Satellite Studies (CIMSS) at 
the University of Wisconsin-Madison has been devel- 
oping prototype methodologies for the assimilation of 
satellite data and products into an analysis/prognosis/ 
display system. The vehicle for the investigations is an 
analysis/prognosis system developed by the Australian 
Numerical Meteorological Research Centre ( ANMRC) 
which has been adapted to the United States region 
and implemented on the Man-computer Interactive 
Data Access System (McIDAS) facility at the space 
Science and Engineering Center. 

Experiments — to date principally done in a research 
mode — have involved techniques for the assimilation 
of soundings from the V1SSR Atmospheric Sounder 
(VAS). cloud and water vapor drift winds and satellite- 
derived surface energy balance descriptors into the 
analyses and prognoses. 

Several authors have experimented with VAS data 
as input to numerical prediction models in both Ob- 
serving System Simulation Experiments (OSSEs) and 
in real case studies. In general, the positive impact of 
VAS from these experiments for the relatively data rich 
continental United States has not been as significant 
as originally had been hoped. The experiments have 
been very helpful, however, both in developing a better 
understanding of the data characteristics and also in 
motivating improved assimilation techniques. It is fair 
to say that progress is being made. 

While comprehensive investigations of the VAS data 
(Jedlovec, 1985: Fuelberg and Meyer, 1986) have in- 
dicated meso-alpha-scale horizontal gradient infor- 
mation content in the soundings, their reduced vertical 


resolution, incomplete spatial coverage and certain 
problematic data characteristics (notably temperature 
biases) have necessitated new approaches to assimila^ 
tion. In general, these approaches are designed to pre- 
serve the information content of the data and model 
prediction fields wffiile circumventing problems inher- 
ent to the observing system. 

One of the formidable limitations of the data source 
is dictated by its incomplete horizontal coverage (no 
soundings in cloudy regions). Gal-Chen et al. (1986) 
demonstrated in a set of OSSEs that the data gaps tvp- 
ical of infrared sounders significantly reduced the im- 
pact of VAS soundings in numerical forecasts and sug- 
gested a geostationary microwave instrument would 
be helpful in resolving this problem. In their experi- 
ments, it was found that multiple-time data insertion 
and nudging of the data into the model were instru- 
mental in achieving forecast impact and reducing 
model shock. 

This study is specifically aimed at another of the 
problematic characteristics of V AS data, the biases in 
the temperature (geopotential) retrievals which are sit- 
uation-dependent and can make their assimilation by 
standard means quite difficult. In these experiments, 
we compare two methods of assimilation. Both are 
variational, but the first is “standard in the sense that 
scalar VAS thicknesses are used as input to the varia- 
tional analysis methodology. In the second method- 
ology. horizontal gradients of thickness derived from 
the VAS soundings are input to the variational analysis 
scheme. 

Cram and Kaplan (1985) have also reported on a 
variational technique for assimilating VAS gradient 
information. The methodology described here differs 
in several respects from their study. Our primary anal- 
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ysis variable is geopotential height, rather than tem- 
perature, to facilitate the later inclusion in the analysis 
of horizontal geopotential gradients inferred from 
measured winds. Also, this variational procedure in- 
cludes vertical structure terms provided from the fore- 
cast model (guess field for analyses); the aim is to pre- 
serve coherent model-generated vertical detail. 

A general problem with the evaluation of satellite 
information is finding verification data with adequate 
space and time resolution to compare. The investiga- 
tion day of 6 March 1982 was one of several AVE/ 
\ AS special network days for whs h a dense space- 
time network of RAOB observations was available in 
the Texas region for verification (Hill and Turner, 
1983). Also, the VAS temperature data at 1445 UTC 
on this day exhibited very strong bias characteristics, 
allowing a meaningful comparison of assimilation 
techniques. 


2. Objective analysis scheme and forecast model 


*-//[■ 


A(A<i> - A<j>) 2 


+B ( d ^- d -^\ + A d -±Jl\ l 

\ch: dxj \<3>' dyj 


\dxdy ( 1 ) 


where 


Fc gross blending functional 
A, B, C reliability weights 

4>, <t> adjusted geopotential and observed 250 mb 
geopotentials (gpm) 

A 0, A 0 adjusted and observed 1000-250 mb geo- 
potential thickness (gpm) 
x, y * longitudinal and latitudinal directions. 


Once this 1000 to 250 mb thickness has been estab- 
lished, individual level geopotential analyses are then 
performed by blending the geopotential gradient anal- 
ysis at a level with the analyzed thickness immediately 
above and below that level. Thus, at any interior level", 
the variational functional is of the form 


a. Analysis methodology 

The three-dimensional objective analysis scheme 
used is an evolution of the scheme developed in the 
Australian Numerical Meteorology Research Centre 
(ANMRC) by Seaman et al. (1977) for limited area 
analysis over the Australian region and adapted to the 
l nited States region. Details and examples of its use 
in this region can be found in Mills and Havden (1983) 
and LeMarshall et al. (1985) and will only be reviewed 
here as it pertains to this experiment. The scheme 
combines the successive correction method (SCM) of 
Cressman (1959) and the variational blending tech- 
niques of Sasaki (1970). 

The analysis system is based upon principal levels 
at 1000 and 250 mb, and the calculus of variations is 
used in an explicitly three-dimensional manner. Esti- 
mates of the geopotential 0 and the geopotential gra- 
dient in each coordinate direction are combined at each 
grid point in the three-dimensional analysis domain. 

The basic observational ingredients combined in the 
blending are observed geopotential thicknesses of each 
layer and horizontal geopotential gradients either from 
geopotential measurements or independently derived 
from wind measurements and a gradient wind law. In 
this study, however, no wind measurement inputs to 
the gradient fields are included. 

A 1000 mb geopotential field is defined by a SCM 
analysis using observations of geopotential and a guess 
field derived hydrostatically from a mean sea level 
pressure analysis and an assumed lapse rate. The 
scheme then uses a three-dimensional variational 
blending of 1000-250 mb thickness and geopotential 
gradient information at 250 mb to produce a bulk 
1000-250 mb atmospheric structure. The variational 
functional for this gross blending is of the form 


Fi~f — A<j> A ) 2 + E(A<p B - A<i> B ) 2 

HS-S'l 


Jd<t> <30 

\a.t s x 


dxdv (2) 


where 


Ft interior blending functional 

D, E, G reliability weights 
A 0 . 4 , A 04 adjusted and observed thickness of layer 

above an interior level (gpm) 

A 0«, A 0 fl adjusted and observed thickness of layer 

below an interior level (gpm). 

This procedure is carried out iteratively in the vertical 
with fixed 1000 and 250 mb geopotentials from the 
preceding gross blend as upper and lower boundary 
conditions. 

Reliability weights (A-G) for the analyzed geopo- 
tential and geopotential gradient fields are determined 
on the basis of data density and latitude and level-vary- 
ing standard deviation parameters. In this way. the 
blended geopotential fields reflect a greater weight to 
the wind data (when used) in low latitudes and to the 
thickness information at higher latitudes. 

The same type of blending is then carried out for 
levels above 250 mb with the exception that onlv the 
2 Mi m level is fixed. The Euler-Lagrange type equations 
which arise from the variational functionals (Equations 
1 and 2 ) are solved through standard over-relaxation 
techniques. 


b. Forecast model 

The Subsynoptic Scale Model (SSM) is a ten-level, 
67 km resolution primitive equation model written for 
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the Australian region (McGregor et al., 1978) and 
adapted for use in the North American region. The 
model in the form for the United States has been 
amended with a more complete boundary layer pa- 
rameterization (Diak et al., 1986), the inclusion of 
Corby differencing (Corby et al., 1 972) to better handle 
the high topographic regions of the United States and 
the replacement of the original Arakawa-Schubert 
convective parameterization with a modified Kuo 
scheme (Kuo, 1965; Hammarstrand, 1977). 

3. VAS data 

The success of the VAS experiment as proposed by 
Suomi et al. (1971) has been the result of close collab- 
oration between scientists and engineers from the Uni- 
versity of Wisconsin, NASA, NOAA Santa Barbara 
Research Center and Westinghouse. The first in a series 
of instruments was put into geostationary orbit aboard 
the GOES-D satellite (Smith et al, 1981). 

The VAS is a radiometer possessing eight visible 
channel detectors and six thermal detectors which sense 


atmospheric radiation in 12 spectral bands whos 
wavelengths lie between 3.9 and 1 5 ^m. A dwell soum 
mode of operation was designed to enable the inter 
pretation of the vertical temperature and moistur 
structure of the atmosphere. Typical space and tirrn 
resolution (in clear areas) are 70 km and 3 h, respec 
lively. Vertical resolution averages about 4 km fo 
temperature and moisture. 

Soundings for 6 March 1982 were prepared with i 
physical retrieval scheme (Smith, 1983) using an LFN 
forecast to provide a first guess. 

The characteristics of VAS soundings for this da> 
have been extensively discussed by Jedlovec (1985) and 
Fuelberg and Meyer (1986) and compared to concur- 
rent RAOB soundings from special network stations. 
Figure la, b shows the mean temperature and geopo- 
tential error structure (satellite vs RAOB) as a function 
of pressure for our investigation time of 1445 UTC. A 
low-level warm bias (negative values) is evident from 
the surface to about 800 mb in the satellite retrievals. 
Above this level a strong cold bias exists to about 600 
mb where a moderate warm bias takes over in tropo- 


MEAN DIFFERENCE (RAOB-SAT) 

6 MARCH 1982 1445GMT 

TEMPERATURE (K) HEIGHT (M) 




Fig. I. (a) Mean error structure of temperature <K) for VAS vs special network 
RAOB for 1445 UTC 6 March 1982: <b ) as in (a) except for geopotential (gpnt). 
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pause levels. Geopotential thickness errors (thickness 
derived from the temperature profiles) may add or 
cancel depending on the vertically integrated temper- 
ature profile. For 1455 UTC we see in Fig. 1 that the 
dominant feature in geopotential is a negative bias 
(positive RAOB-satellite values) from 800 to almost 
300 mb. 

4. Synoptics — 6 March 1982 

Figure 2a shows a conventional 500 mb RAOB 
height analysis (also analysis and forecast model do- 
main) for 1200 UTC 6 March; and Figs. 2b and c in- 
dicate the corresponding surface and 700 mb dewpoint 
analvses, respectively. Conditions for this day have been 
discussed by Jedlovec (1985). The case is marked by a 
developing frontal wave in the southeast United States 
with a significant thickness and dewpoint gradient back 
through the Texas region. The system was the cause 
of convective activity in the southeast later in the day. 
Behind this frontal system, a shallow high pressure sys- 
tem was present around western Oklahoma. To the 
north, a cold front from a low pressure system in Can- 
ada was pushing south through the upper plains. 


5. The experiments 

Two analyses using VAS data were made at the 1445 
UTC VAS observation time on 6 March. The first 
analysis was made incorporating only scalar VAS 
thicknesses as input into the variational analysis 
scheme. No horizontal geopotential gradient infor- 
mation, either from VAS measurements or derived 
from wind information, was used in this analysis. 1 he 
second analysis for this time was made using only hor- 
izontal geopotential gradient information from the 
VAS as input to the variational functionals (Equations 
1 and 2) and no scalar geopotential information. Dew- 
point temperatures for these two analyses were iden- 
tical. both being standard SCM analyses from the scalar 
VAS information. An analysis was also done at nuu 
UTC 6 March using the special network RAOB 
soundings, to be used as verification for the 1445 UTC 
VAS analyses. The first-guess fields for the two V AS 
analvses and the special network analysis were PfO'-ided 
by an SSM forecast to 1 500 UTC made from a RAOB 
analysis at 1 200 UTC (Fig. 2). f 

A control forecast of 12 h duration was made from 
the 1200 UTC RAOB analysis. Two forecasts were 
made from the 1445 VAS standard and gradient anal- 
vses to the same time (0000 UTC 7 March). A third 
forecast, 1445 to 0000 UTC, was also made using the 
gradient geopotential analysis with the first-guess ( 1 ^0 
UTC control forecast) moisture field to assess the rel- 
ative roles of mass versus moisture information on 
precipitation simulations. Lastly, a RAOB analysis was 
also done at 0000 UTC to be used as verification tor 
these forecasts. 



S. Results 

a. I 'AS analyses 

Figure 3a-d shows the VAS standard and gradient 
analvses at 1445 UTC, the 500 mb special network 
RAOB analysis at 1 500 UTC, and the guess field ( 1 500 
UTC control forecast) for these analyses. Data points 
on Fies. 3a and b are VAS retrieval locations, while 
those on Fig. 3c are special network sounding locations. 
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Fig. 3. (a^ The VAS standard 500 mb analysis and VAS sounding location 1445 UTC 6 March 1982; (b) 
VAS gradient 500 mb analysis; (c) 500 mb special network analysis 1500 UTC 6 March 1982 and special 
network sounding locations; (d) 500 mb 1500 UTC control simulation (first-guess field). 


Some of the potential problems in incorporating 
VAS soundings into synoptic analyses have been dis- 
cussed by Cram and Kaplan (1985) and others. These 
include bias errors and potential discontinuities be- 
tween data and no-data regions. Both of these problems 
are seen in the standard VAS analysis at 1445 UTC. 
In the region of the special network data coverage, we 
see that the VAS standard analysis (Fig. 3a) has 500 
mb heights which are biased as much as 30 m low r 
compared with the special network analysis (Fig. 3c). 
Some evidence of discontinuity (problems of the anal- 
ysis model in reconciling the data with the guess field) 
is evident in the Colorado region. The problem is ex- 
acerbated by the fact that the first-guess field (Fig. 3d) 
shows some positive geopotential bias at 500 mb com- 
pared to the special network analysis. A stronger dis- 
continuity is evident in the southeastern United States, 
where the combination of biased VAS data and poor 
data distribution presents a difficult analysis problem. 

Figure 3b is the VAS gradient analysis at 1445 UTC. 
The discontinuities evident in the VAS standard anal- 
ysis have been eliminated Here, perhaps at the price of 
somewhat oversmoothing this 500 mb height field. The 
gradient analysis has modified the guess field in the 
correct direction in the special network region, reducing 
the positive bias. Gradients appear fairly well repre- 
sented although the gradient analysis has produced a 
flattening of gradients in central Texas and a smoothing 
of some of the trough features in the Nebraska-Okia- 


horaa region. Jedlovec (1985) has noted a general re- 
duction in horizontal gradients from the VAS com- 
pared to special network RAOB data for this day, w hich 
is a partial cause of the problem. Also contributing is 
the inherent smoothing in calculating gradients over 
finite grid element lengths. A short wave feature has 
been created in the southeast United States in the gra- 
dient analysis w hich is outside of the special network 
domain and cannot be verified. 

Figure 4a, b shows differences, respectively, between 
standard and gradient analyses and the special network 
analysis for the special network verification domain. 
As shown in this figure for both standard and gradient 
analyses, the maximum errors in the verifiable region 
w'ere in the east, where there was a combination of a 
large VAS data gap and a maximum in the guess field 
error. 

Table l shows biases, standard errors and Si skill 
scores for the standard and gradient analyses versus 
the special network analysis over the special network 
domain. As seen in this table and visually in Fig. 4, 
the gradient approach to VAS data assimilation has 
greatly reduced the errors produced by strongly biased 
information. 

b. WAS forecasts 

Forecasts of 500 mb geopotential valid at 0000 UTC 
7 March are compared in Fig. 5a, d. Figures 5a and b 
are forecasts from the VAS standard and VAS gradient 
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Fig. 4. (a) Difference between VAS standard analysis and special 
network analysis at 500 mb; (b) difference between VAS gradient 
analysis and special network analysis at 500 mb. 


analyses, respectively, done at 1445 UTC, and Fig. 5c 
is the control forecast of 12 h duration made from the 
synoptic analysis at 1200 UTC 6 March. Figure 5d is 
the verifying svnoptic analysis for these forecasts at 
0000 UTC 7 March. 

Many of the differences evident between gradient 
and standard analyses have been smoothed out in the 
subsequent forecast process and only more subtle dif- 
ferences are evident. The wavelength of the trough has 
been shortened by both VAS forecasts compared to the 
control and verification. Both have the trough some- 
what too deep in the mid-Texas region, with the stan- 
dard forecast slightly worse, probably due to the neg- 
ative geopotential biases which existed in its initial 
analysis. Geopotentials in the southern Illinois and 
Missouri region appear to be somewhat better in the 
gradient versus standard forecast. 

The control forecast is generally good, but exhibits 
somewhat weaker gradients than the verification along 
the trough axis and some positive bias in the northern 
Texas region. Intuitively, it is difficult to call any one 
of these forecasts superior to the others. 


Biases, standard errors and S 1 skill scores for these 
three forecasts versus the 0000 UTC synoptic analysis 
are shown in Table 2 for a region encompassing the 
VAS data points which are shown in Fig. 3a. These 
statistics bear out the intuitive evaluation of a mixed 
set of results from these forecasts. For clarity of pre- 
sentation, the gradient statistics have been marked in 
Table 2 where they have outperformed the standard 
and/or control forecasts. As shown, the gradient fore- 
cast dearly outperforms the standard forecast, with the 
most impact being in S 1 skill scores and standard error 
statistics. In the gradient versus control forecast, im- 
provements are noted in biases at many levels, while 
standard error comparisons between the two are more 
mixed. The S 1 skill scores between gradient and control 
forecasts are always within about one percentage point 
of each other, however, at no level is the gradient fore- 
cast able to outperform the control in this comparison 
with a verification analysis made using synoptic-scale 
data. 

Other investigators (Gal-Chen et al., 1986; Cram and 
Kaplan, 1985; Mills and Hayden, 1983; and others) 
have similarly found small impact of subsvnoptic scale 
satellite soundings on primary forecast fields for the 
continental United States. Mills and Hayden (1983), 
however, found significant impact on forecasts of pre- 
cipitation, dewpoint and vertical motion and it thus 
appears worthwhile to examine these fields and their 
interrelationships also. 

Precipitation patterns for the three forecasts for land 
areas are shown along with verification in Fig. 6a-d, 
and it can be seen that indeed some significant differ- 
ences exist. Precipitation in the SSM is produced by a 
combination of large-scale supersaturation removal 
and/or convection parameterized with a Kuo-type 
scheme (Kuo, 1965). In prior sensitivity experiments 
(Diak et al., 1986), changes in convective precipitation 
for this model between similar forecasts were shown 
to be most dependent on changes in surface moisture 
convergence, which is the forcing term for the Kuo 
parameterization. In the 6 March experiments de- 
scribed here, while there were some areas of large rel- 

Table l. Biases, standard errors, and SI skill scores for 1445 
UTC standard and gradient analyses vs 1500 UTC special network 
analysis. 

Standard analysis Gradient analysis 


Level 

(mb> 

Bias 

(m) 

SE 

(m> 

SI 

Bias 

(m) 

SE 

(m) 

SI 

1000 

2.21 

7.65 

31.18 

2.24 

7.59 

30.41 

850 

6.25 

23.00 

65.99 

2.32 

9.17 

33.32 

700 

-*1.38 

20.14 

40.99 

■ 0.57 

10.36 

25.57 

500 

-13.69 

19.82 

26.15 

-1.59 

10.35 

12.96 

400 

-23.84 

20.36 

24.88 

-2.17 

12.38 

12.71 

300 

-32.66 

22.62 

24.82 

-2.14 

14.29 

11 44 

250 

-26.63 

26.78 

26.75 

-0.63 

14.81 

10.32 

100 

-29.95 

29.52 

28.P 

-2. S3 

14.17 

11.69 

150 

-49.61 

30.14 

33.80 

-4.06 

14.09 

1 3.65 

100 

-60.48 

31.12 

46.33 

-3.99 

13.29 

17.89 
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ative changes of convective precipitation between fore- 
cas s, it was large-scale precipitation which dominated 
both regional and overall model totals and produced 

diR>rPn^ S St ^ W !l m F,g - 6 ' Accordin gJv. we look to 
differences in the dewpoint and vertical motion fields 

as being responsible for these precipitation changes 
Vertical motion fields at 1800 UTC (* at a = 0.49 
approximately 500 mb for no terrain) are shown in 
big. 7 for the two satellite forecasts and the control 

f ° r thC forec ^ts, the 

500 UTC « uess fie!d (control forecast) for this analysis 


and the 1500 UTC special network dewpoint analysis 
are shown m Fig. 8a-c. ' 

Precipitation patterns from the control forecast (Fig 
6c) are reasonable, but several deficiencies are evident' 
^^ofUke Erie, large amounts of precipitation are 
predicted where verification shows a smaller area of 
lesser accumulation. Precipitation along the southeast 
front is not as well defined as the verification indicates 
and maxima are reduced and somewhat misplaced 

Ihe standard forecast shows a minor reduction in 
tne northeast maximum predicted in the control fore- 


_]ABLE2_ Bi a ses. standard errors, and SI skill scores for standard. 
Standard fcst 


1000 

S50 

5i*0 

400 

300 

250 

200 

150 

100 


gradient and control forecasts vs verifying 0000 UTC 


synoptic analysis. 


Gradient fcst 


29.67 

19.98 
5. IS 

-6.54 

-17.96 

-7.49 

3.42 

0.75 

13.99 


Control fcst 



20.03 
16.76 
16.11 
16.63 
24.27 

34.03 
33.29 
31.45 
28.42 
27.60 


79.29 
66.32 
37.57 
21.38 
23.06 
21.69 
2 1.15 
22.01 
24.37 
30.11 


1 7.38** 
2I.6T* 
13 . 01 ** 
-0.1 6** 
-8.85 
-18.96 
-7.78 
5.13* 
5.68* 
23.73 


^Outperformed standard forecast. 
Outperformed control forecast. 


17.89* 
16 . 10 * 
14.55** 
18.06* 
23.65** 
32.7-* 
30.96** 
30.26* 
25.52* 
20.85* 


68.26* 

58."6* 

31.54* 

20.91* 

20 . 88 * 

19.80* 

18.92* 

19.39* 

20.67* 

22.93* 


17.56 

22.06 

17.18 

11.27 

2.63 

- 8.81 

0.98 

11.92 

10.13 

25.98 


15.57 

15.23 

16.27 

24.52 

34.33 

33.62 
31.79 
24.08 

19.62 
9.99 


63.28 

54.39 

31.14 

19.32 
19.46 

18.56 
17.81 
18.50 

19.33 

21.56 
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, FlG - 6 ; < a) The \ VAS s^ndard forecast of precipitation: (b) as in 6a except gradrent 
totajT ’ C ^ m 63 CXCePt C ° ntr01 forecast; (d) l200 -° 000 UTC Service A precipitation 


cast, but the frontal precipitation in the southeast is 
grossly overpredicted and the axis has been unrealisti- 
cally elongated due to a spurious maximum in Loui- 
siana through east Texas. As can be seen in Fig. 7a, 
the precipitation in this area exactly coincides with a 
large predicted area of intense vertical motion. Recall 
that in the analysis used to initialize this forecast, this 
was an area of discontinuity caused by data biases and 
guess-field error. It is likely that mass adjustments in 
this presumably sensitive cold frontal region have ini- 
tiated the vertical motion and precipitation. Evidence 
ot these adjustments was seen in large surface pressure 
oscillations in the region in the first several hours of 
this forecast. Similar to Cram and Kaplan ( 1 985) find- 
ings, these oscillations were significantly reduced in the 
gradient assimilation. 

The gradient forecast of precipitation (Fig. 6b) is, in 
general, acceptable with the frontal precipitation 
sharply defined and bearing a close resemblance to the 


verification in this region. The vertical motion patterns 
and precipitation do not show the problem associated 
with the standard forecast and also appear more well 
defined than the control simulation. In the northeast, 
however, the overprediction of precipitation has been 
only slightly reduced from the control run. 

To better clarify the roles of moisture and vertical 
motion in the changes in precipitation between fore- 
casts. one additional 1445-0000 UTC simulation was 
made (gradient- 1 experiment) using the gradient geo- 
potential analysis at 1445 UTC, but retaining the 
moisture fields from the first guess ( 1 500 UTC control 
forecast. Fig. 8b). Thus, this forecast has new satellite 
mass, but no moisture information. 

In the frontal region, the precipitation and vertical 
motion patterns from this new forecast very closely 
resembled that of the gradient experiment. It is evident 
from a comparison of the 1445 UTC satellite moisture 
analysis with the guess field and the special network 


nod jo 



Fl<i. 7. (a) vertical motion (<> at a - 0.49. approximately 500 mb) 1800 UTC from VAS standard forecast (sec -1 X 10‘, negative upward); 

(b) as in 7a except gradient forecast; (c) as in 7a except control forecast. 
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analysis moisture (Fig. 8), that the satellite information 
did improve moisture gradients in this region. The gra- 
dient- 1 experiment suggests, however, that it was the 
mass information from VAS and the subsequent ver- 
tical motion prediction which were responsible for the 
well-defined precipitation patterns in the original gra- 
dient forecast. 

Although VAS soundings in the northeast were 
sparse and poorly distributed, evidence points here to 
VAS dewpoint information as the cause for the small 
reduction in the precipitation maximum in the stan- 
dard and gradient forecasts versus the control. As 
shown in Fig. 8, the dewpoint values as well as hori- 
zontal gradients in the region have been somewhat re- 
duced in the VAS analysis versus the guess field. The 
0000 UTC synoptic analysis showed a similar weak- 
ening of dewpoint gradients. While the gradient forecast 
had upward vertical motion which was as strong as the 
control in this area (Fig. 7), there was still a reduction 
in precipitation in this forecast, suggesting predomi- 
nance of the dewpoint information. Additionally, the 
gradient- 1 forecast, without any satellite moisture in- 
formation, displayed precipitation in the region which 
was very similar to the control simulation, adding ev- 
idence that moisture information was responsible in 
this case for the precipitation changes between exper- 
iments. 

7. Conclusions 

A comparison has been made between analyses and 
forecasts which incorporate VAS geopotential infor- 
mation by scalar and gradient variational assimilation 
methodologies. For this case of 6 March 1982, a day 
on which the VAS data exhibited strong biases in geo- 
potential, the gradient methodology clearly outper- 
formed the standard in error statistics for the 1445 UTC 
comparison of analyses to special network RAOB ob- 
servations. The gradient forecast outperformed the 
standard scalar methods in most statistical comparisons 
and was about the same as a control simulation made 
w ithout satellite data for overall quality of the geopo- 
tential forecast. Effects were noted in the gradient fore- 
cast versus both the standard and control forecasts in 
vertical motion patterns and the prediction of precip- 
itation. 

The case is especially encouraging because it is a 
“worst case scenario” for VAS data — namely, a case 
in which ther : were very significant data problems, but 
in which appropriate assimilation methodologies re- 
sulted in an analysis/forecast which was, at worst, not 
degraded from a control comparison and in some re- 
spects was improved. Further work will be aimed at 
concurrently optimizing the gradient assimilation 
methodology and the VAS sampling strategy to extract 
the most information from the observing and assimi- 
lation systems. 

Results from this study suggest that a potentially 
powerful observing system may be a combined RAOB 
and satellite system, with the RAOB observations used 


to establish synoptic-scale scalar values of observables 
and satellites defining the mesoscale gradients. 
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1. 


INTRODUCTION 


System Interactlva Data Access 

in the 107 ( 1 ' ^ ’ 3S itS name developed 

the 1970 s as primarily a satellite data ? 

access and display system. Kotable contribu- 
ons w«re cade in areas of cloud and water 
vapor drift winds, nowcasting products, satel- 

““ h g$ ’ Car ! h radlation budgets, plane- 

tarj atmospheres and many others. The potential 

.e «_.r ginal system, however, was limited by 
its computer architecture; that is several 
smaller computers sharing a data base with each 
driving several display terminals. Recentlythe 
research emphasis has been expanded with a 

f " Tr\ ter strate £T to ore mainframe CPU 

(IE.. 4381) driving all McIDAS terminals 

a 1 2 j ™ 5 a y ailab IHtv of the mainframe has 
allowed expan sicn into applying some of the 
research products for synoptic and mesoscale 

eoai y is% an a f0reCasCln 8- A Primary research 
t0 determine the potential benefit of 
satellite data and derived products for such 
applications. To date, research efforts have 
ncentroted on techniques for assimilation of 
If) 1 soundings , cloud drift and water vapor 

b-lfnr W1 H S ** saCellite- derived surface energy 
balance descriptors in numerical analvses and 

prognosis. Concurrent with this research effort 
5 the development of a prototype M 4-d" 
assimilation scheme to examine’ the feasibility 
cf regular updating of a forecast model with ’ 

g°a 1 ^ mT 't T" aS ^°P cic information; the 
E^al being the best possible continuous 

f r !f r ,^ 1C T' cf at=os Pheric state combining data 
r "‘“ aU observing systens. 

DATA AND PRODUCTS 

* * * VAS Soundings 

Sou-der I!:? aUCCeSS ° f riSSR Atmospheric 

a 1 C’q-m f e ^P e niment as proposed by Suomi et 
Co ii.- ' has been th<? result of close 

Deration between scientists and engineers 


slnla ° f Wiscpoa ln. NASA. NOAA 

fi£t E* Research Cen£er and ^Inghouse. 
1 / 6erlCS ° f lnstru »«nts was put 

-see iw° St rr l0nary orbit aboard the COES-D 
satellite (Smith et al. p I9fll), 

, , * ' AS * s a radiometer possessing 

e g t visible channel detectors and six thermal 

iHjectral V'T a£,DOS P herlc radiation in 
V, f Whose wavel «°gths lie between 

tion (in rl 1 ”’ Typl f al ®P ace and time resolu- 
tion (in clear areas) are 50 km and three 

hours, respectively. The soundings are pro- 

fll't ^ r 1 tlBe “ C the ^ 0o P erat ive Institute 
thl ” e . te ° ro l°8 lcal Satellite Studies (C1MSS) at 
the University of Wisconsin using a full 
physical iterative solution of the radiative 
transfer equation. The general qualities of 

« m ri ar (!5L)° Un J 1, J 8 r, are d€scribed in Velden 
et al. (1984) and Jedlovec (1984). 


2.2 


Cloud and Water Vapor Drift Wind; 


„ f loud and water va P°r drift winds are 
also produced at the CIMSS on the McIDAS 
faculty using the technique of Mosher (1978). 
This data base is generated from GOES 
half-hourly visible and infrared imagerv and 
hourly water vapor imagery. Wilson and 
Houghtcn (1979) have estimated the standard 
C < rr j r ° t J the dlf (ertnce between cloud drift 
winds and ravinsor.de reports at about 4.7 ms' 1 

Sta£ i stlcs water vapor drift 
w.r.ds (..osher and Stewart, 1981 ) show a 
standard error of about 8 cs~*. 

yb e height assignment of cloud drift 
winds h«vily on the bi-spectral method 

^ Heigh£ assl ^ ent water 
vapor w.nds (mid-troposphere level) relies on 
the use of water vapor brightness temperatures 
and nearby environmental soundings. 


2.3 


.S urface Ener gy Balance Descriptors 


The description of surface energv 
budget parameters relies on the time change of 
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the VA5 radiances or soundings from the lower 
atmosphere being a descriptor of the sensible 
heating at the surface (l.e., balance of 
sensible/latent heat). Thi6 information, 
calibrated through use of a surface-boundary 
layer model, is used to derive surface energy 
balance descriptors for prognosis models. Fig. 
1 shows such a calculation and represents the 
non-dimensional ratio of sensible heat to net 
radiation at the surface both calculated from 



ratio of surface sensible heating to net 
radiation July 20, 1981.(2 X 100) 


3. THE ANALYSIS SCHEME AND PROGNOSIS MODEL 

The principal characteristics of the 
analysis scheme and forecast model used in this 
study are shown in Table 1. 

3 . 1 Analysis Scheme 

The analysis scheme uses a 51x61 point 
grid with a 67 km spacing, on a Lambert 
Conformal projection covering most of the 
continental U.S. The scheme is an adaptation to 
the North American region of the scheme written 
by R. S. Seaman at the Australian Numerical 
Meteorology Research Centre and is described in 
Mills and Hayden (1983). The analysis scheme 
uses Che successive correction technique of 
Cressman (1959) to generate grid point fields 
frcm observations. Subsequent to this the 
variational blending techniques of Sasaki (1958, 
197C) are used to combine fields of geopotential 
height, thickness and gradient (wind), weighted 
by their reliability, into an "optimal" 
aralvsis. The scheme uses an initial 
three-dimensional blend of 10CO mb and 250 mb 
height, gradient and thickness data to produce a 
bulk atmospheric structure. 

Individual level geopotential analysis 
nre then performed by blending the geepctent ial 
gradient analyses at a level with the analyzed 
tbickr.ess immediately above or below* that level. 
This scheme is veil-suited to the analyses of 
satellite temperature profiles ar.d w*ater vapor 
and cloud drift v-inds because gerpctencial 
thickness is a primary analysis variable, 
tacilitating the assimilation of satellite 
temrerature profiles into the objective 
analysis, while the variational blending ensures 
a consistent introduction of gradient (wind) 


information into the analysis. The blended 
geopotential fields and the gradient/geostrophic 
relationship are used to produce a wind field 
from which a stream function field is calculated 
by a variational blending procedure (Seaman et 
al., 1977). The forecast model commences its 
Integration from this non-divergent field. 

3.2 Prognosis Model 

The forecast model is a ten level, 67 km 
resolution primitive equation model, originally 
written for the Australian Region (McGregor et 
al., 1978) and then adapted for use in the North 
American Region (Mills et al. v 1981). The 
Australian Region version of the model is used 
operationally by the Australian Bureau of 
Meteorology. The North American version of this 
model is different in several respects. It uses 
the finite differencing scheme of Corby et al. 
(1972) to minimize truncation errors over steep 
topography. It uses a Kuo type convective 
parameterization (Kuo, 1965, 1974) which 
replaces the Arakawa-Schubert scheme described 
in McGregor et al. (1978) and it also has a 
comprehensive boundary layer scheme (Diak et 
al., 1986) using a similarity-theory surface 
layer coupled with vertical diffusion of heat, 
moisture and momentum at the three lover model 
levels, surface shortwave and longwave radiation 
modified by cloud and a surface energy balance 
evaluation. 

Table 1 

Principle Features of the Analysis Scheme 

and Prognosis Model 

Analysis Scheme 

Combination of successive correction method and 
variational blending in three dimensions. 

Ten pressure levels, p - 1000, 850, 700, 500, 

400, 300, 250, 200, 150 and 100 mb. 

Horizontal resolution: 67 km 

Fields analyzed or derived at each pressure 
level : 

- geopotential height 

- temperature 

- dewpoint 

- wind components 

- stream function 


Prognosis Model 

Primitive equations model in sigma coordinates. 

Ter. vertical levels at * 0.C9, 0.19, 0.29,..0. CQ 

Horizontal resolution: 67 kn 

Staggered horizontal (Arakava "C") grid 

Semi-irp licit tire differencing 

Sur: ace/boundary laver parareterizat ion 

Sroad scale prec ip i tat ion 

Kuc-type convective pa rare r e r i za t ion 

Updated boundarv conditions 
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4 . INFORMATION POTENTIAL 

4.1 VAS Soundings 

Fig. 2 (a-c) depicts a aeries of 500 mb 
geopotentlals for the AVE/VAS day of Hatch 4 
1982 when high density space and time radiosonde 
data was available ir. the Texas region for con- 
firmation of analyses and prognoses. rl &* ‘ . 

shows an analysis (with AVE sounding lotions) 

Ba de from this radiosonde daca atl5 CKT. Fig. 

2b Is an analysis made for 1430 GMT readjusting 
the guess field (Fig. 2c) to match horizontal 
height gradient Information from the VAS sound- 
ing! In a version of the variational analysis 
methodology previously outlined 
points are also depicted on Fig. 2b). The gues 
field (Fig. 3c) for both analyses Is a forecast 
to 15 GMT made from a conventional network 



Fig. 2 (aY! Special network 500 mb analysis 



Fig. 2 (b). Same as (a) only VAS gradient 



Comparing these figures, we see that the 
VAS gradient analysis has corrected sos« of the 
guess field's high geopotentlal bias within the 
data domain and brought the fielda Into b «tter 
agreement with the special network radiosonde 
snslysls, perhsps st the price of some over- 
smoothing of the field. While the VAS soundings 
themselves for this time were significantly 
biased, the gradient assimilation methodology 
has greatly reduced these problems for analysis 

of the data. _ . 

A forecast frcu this VAS analysis 

(1430-00 GMT) demonstrated similar skill as a 
control forecast (1200-00 CKT) in geopotentlals. 
Improvements in precipitation patterns were 
noted in the VAS forecast versus the control. 

4.2 Cloud and Water Vap or Drift Winds 

Fig 3 (e-d) shows an analysis situation 
particularly Impacted by the Input of cloud and 
water vapor drift winds for s case study day 
April 26, 1982 (Le Marshall et si.. 1986). 
Inspection of satellite imagery for this day 
showed the presence of a small cut-off ° 

500 mb which developed In the southern Illln 
region between synoptic stations and symopt c 
times at about 1600 GMT and gradually weakened 
during the course of the day. The VAS *°"“ di ^ s 
alone, done only in clear areas, were unable to 
depict this feature and it did not develop in 
control forecast (1200-00 CKT) made from 
' synoptic data. Using cloud and water vapor 
drift winds, however, in analyses we were able 
to capture this feature and, in part, its time 

evolution^ 26> Flg . 3 ( a -d) shows VAS 

gradient plus water vapor and cloud drift * i !L_ 
analyses Lde for 1600, 1700, 1900 and 2000 CMT 
on this day with the satellite winds also shown. 
The influence of the wind data is evident an 
principle it can be seen that the satellite 
sounding and wind data have complementary 
coverage with the former available only in clear 
areas and the latter mostly in cloudy regions. 

5. THE 4-DIMENSIONAL ASSIMILATION STSTFM 

One of the assets of the McIDAS system 
is the efficient manner in which conventional 
and satellite data products are routine. y 
ingested and processed. One of the greate 
changes which has come about in recent years is 
the huge increase in volume of satellite 
products as the spatial and tempera sampling 
frequency of these instruments has increased, 
and interactive processing techniques have 
refined. The challenge is to provide the be t 
possible four-dimensional specification o. e 
state cf the atmosphere, utilizing the 
information content of these differing 
information sources, and allowing fer^t e 
particular strengths and weaknesses o. eac . 
kat is. to allow the different data types to 
complement each other. This problem will be 
further exacerbated in the near future when an 
operational ground-based wind profiling network 
is implemented over the U.S.. which will provide 
sneci'ication of the vertical wind profile with 
a high degree of accuracy, and with a hitherto 
unavailable sampling frequency. 

A four-dimensional data assimila 
scheme has been implemented on the McIDAS 
svstem. with the aim of (1) providing gridded 


319 



ORIGINAL RAGS fS 

OF POOR QUALITY 




Fig. 3 (a-d). Analyses performed for 1600, 
1700, 1900 and 2000 CMT April 26, 1982 with VAS 
gradient and cloud and water vapor drift wind 
input . 


fleld6 of primary atomapheric variables, 
utilizing-all data available to the McIDAS 
system, for diagnostic studies, particularly 
during special effort field experiments such as 
GALE and STORM; (2) to provide guess fields for 
satellite retrieval processing during the 
sssiflilatlon; (3) to provide a vehicle for 
further investigation of the problem of 
assimilation of data from nixed sources at 
mesoscale resolution; and (4) to provide initial 
state specification for NVP data impact studies 
of significant veather events. 

To achieve these alms, the systen has 
been designed to have a variable grid spacing in 
the neso-a scale (i.e., 50-150 km) , a time 
resolution (that is, analysis update frequency) 
of 3-12 hours, to have a variable domain, and to 
have the ability to nerge data from all sources 
at CIKSS. The chief components of the 
systen are the three-dimensional variational 
analysis scheme described in Mills and Kayden 
(1983), but modified to Include data-density 
dependent influence radii, an elliptical 
influence region for the moisture analysis, and 
an analysis of wind component deviations from 
the background field. The forecast model used 
to link the analyses in time i6 the new ARPE 
model (Leslie et si., 1985). The model is 
nested In a broad-scale forecast, has fully 
parameterized horizontal and vertical dimensions 
and includes a full suite of the eo-called 
"physical parameterlzations,” and a vertical 
mode initialization scheme. Analysis increments 
are added to the model fields at each analysis 
time. 

.Care has been taken to keep the 
individual components of the system as modular 
as possible, so that they may be replaced by 
test modules for assimilation technique 
development, and a number of options are 
available. For example, the particular mix of 
data required is user specified from that 
available in McIDAS; the use of divergent or 
non-divergent winds for model initiation, the 
choice of wind law used to balance the mass and 
wind data, the selection of data rejection 
tolerances, maximum and minimum influence radii 
for each data type, and a vide range of 
prognosis model options are user specified. 

Output products from the system include 
fields of mean sea level and surface pressure, 
surface temperature and dewpoint, and multi- 
level fields of divergent and non-divergent wird 
components, temperature, dewpoint, geopotential 
height, mixing ratio and vertical motion. In 
addition, a vide variety of diagnostic 
quantities can be derived from these fields 
using McIDAS applications software. All 
products are available on the McIDAS systen. 

So far, the scheme has been tested in an 
assimilation code at three hour update freauercv 
using the VAS temperature and moisture sounding 
and gradient w’ind data, the special network 
radiosonde data end the surface data free the 
March 6, 1982 situation, and at a 12-hour update 
frequency during a nine-day period of October 
1985 using VAS temperature and moisture profiles 
and gradient winds, and cloud drift and water 
vaper winds prepared for the KMC, together with 
the radiosonde network data. In addition, a 
number of individual case studies taken from the 
19S5 hurricane seascr. have been analyzed, with 
the domain being centered over the particular 
tropical storm on each day, and the resolution 
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f'i.e., grid spacing) selected to complement the 
particular data distribution of the case. This 
domain and resolution selection is interactive. 
Figures 4 and 5 show examples of analyses 
produced in real time using the system. The 
differing domains were selected on the basis of 
the data distribution and the expected external 
forcing. Figure 4 shows the 200 mb flow pattern 
over Hurricane Elena, where the upper level 
outflow is well-defined by the cloud drift 


winds. Figure 5 shows the low level flow around 
Hurricane Bob, again with the cloud-drift winds 
shown. 



showing the outflow over Hurricane Elena. The 
cloud drift data near that level are overlain. 



cloud drift wind data shown over the Gulf of 
Mexico complemented the analysis of the low 
level vortex of Hurricane Bob. 
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ABSTRACT 

Four simulations are made with a numerical mesoscale model to 
determine the effects of suppressing the initial mesoscale infor- 
mation in the moisture and wind fields on the precipitation fore- 
casts. The second half of a 12-hour control forecast represents 
the ideal reference forecast. This reference forecast can be 
thought of as being "initialized" with fields from the six-hour 
control forecast that contain both mesoscale and synoptic scale 
features. The first experiment is initialized with the same six- 
hour control forecast, but the moisture fields are smoothed to 
eliminate the mesoscale component. The horizontal wind fields 
are replaced with nondivergent winds in the second experiment, 
but all the other variables come from the six-hour control. This 
experiment therefore contains only a topographically-induced ver- 
tical motion field initially. The third experiment contains both 
smoothed moisture fields and nondivergent winds initially and is 
considered the experiment with the most degraded conditions. 
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All of Che experiments' six-hour simulations produce less 
precipitation than the control forecast. Precipitation rates lag 
behind the control forecast, and areal coverage is generally less. 
The third experiment has the worst performance compared to the 
control, which is expected based on its initial conditions. The 
smoothed moisture experiment produces less precipitation than the 
nondivergent experiment in areas over the ocean where moisture 
magnitudes are reduced by the smoothing process. Over a land 
area, however, the smoothed moisture experiment forecasts more 
precipitation than the nondivergent experiment, although moisture 
magnitudes are reduced in that area also. Generally, the differ- 
ences between the first two experiments' forecasts are not large 
enough to conclude whether either mesoscale moisture or initial 
vertical motion is more essential to precipitation forecasting 


with a mesoscale model. 
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A case study evaluation of satellite-derived rainfall 
estimates and their application to numerical model 
precipitation forecast verification 

by 

Glenn Alan Field 

• • Professor David D. Houghton 

under the supervision of Professor 


abstract 


Sat ell i tender ived precipitation estimates are computed 
and then evaluated using a dense network of cooperat 
observer rain gauge reports as the verification. The 
feasibility of using these satellite rainfall estimates to 
evaluate numerical model precipitation forecasts is 
investigated. The correspondence between the numerical 
model forecast and the observations also is assessed. 

The satellite rainfall estimates are produced every 

• j cf artina 1200 GMT/ July 
half-hour for the 24-hour period start g 

20 1981. They are computed using the operatio 

Self ield-Oliver Convective Rainfall Estimation Technigue on 

the University of Wisconsin's Man-Computer Interactive Data 
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Access System (McIDAS). (Suomi et. al., 1983) A severe 
weather outbreak occurred over parts of the southern 
Midwest during this period and significant rainfall amounts 
were observed. More than 300 cooperative observer rain 
gauge observations made during the same time period as the 
estimates are compiled. The McIDAS analysis procedure 
provides estimate values assigned to grid points spaced 22 
km apart. The rainfall observations, however, are at 
irregularly located positions. In order to be able to 
objectively evaluate the estimates, the observations are 
interpolated to the same grid points as the estimates using 
a minimum of smoothing. Difference fields then are 
evaluated . 

The numerical model evaluated is an Australian 
mesoscale model referred to as the Subsynoptic Scale Model 
( SSM ) . Its 24-hour precipitation forecast is examined for 
the same time period as the satellite estimates and 
ground-based observations. The horizontal resolution 
(134 km) and map projection of the SSM are much different 
than for the estimates and observations. A regridding and 
interpolation scheme is employed, which allows the SSM model 
to be objectively evaluated on a common grid with the 
estimates and observations. 


i i i 

The results show that the satellite estimates compare 
very favorably with the observations, especially with 
regard to location of rainfall maxima. It is shown that 
the orientation of the maxima and minima axes in the 
contoured estimate field is in good agreement with the 
observations and radar reports. As would be expected, this 
agreement improves with higher amounts of smoothing. There 
are many apparent overestimates, for which several 
plausible explanations are given. Some displacement errors 
are observed and it is shown how small location errors can 
lead to large errors in a gridded difference field. 

By using satellite estimates as part of the SSM model 
verification, this study suggests a new application for the 
use of the Scof ield-Oli ver technique. Unfortunately, the 
SSM model fails to accurately predict convective precipi- 
tation in this case study. Its forecast precipitation area 
is too far to the north and the amounts are much too 
small. Nevertheless, the feasibility of using satellite 
estimates to verify the model is demonstrated. It is shown 
that the potential exists for operational numerical 
(mesoscale) modeling to benefit by having such satellite 
verification information for precipitation which can be 
produced in near real-time. 
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1. INTRODUCTION 

the Man- computer Interactive Data Access 
Svsten (McIDAS) (Suomi et al. , 1983) 
extensive capability to ingest c ° nv ' n ^°" 3 in 
surface and upper air meteorological data, in 
add •' tion to satellite-derived data which can be 
processed and edited in the McIDAS environment. 

These satellite products rou ^" ely 4™ S *“ h 
temperature and moisture profiles from both 
orbiting (Smith et al.. 1979) and Reostationary 
(Smith, 1983) satellites, as well as cloud drift 
and water vapor winds derived from geostationary 
satellite imagery (Mosher. 1979, Stewart et^al. . 
1965). and gradient winds estimated from VA5 
temperature soundings (Hayden, 

While large amounts of application 
software exists to analyze and process these 
data for specific case studies, including 
advanced objective analysis and numerical 
weather prediction models (Mills and 

1983, LeMarshall et al., 1985), no forma 

mechanism has hitherto existed 

routine assimilation of these mixed data types 
with the aim of producing four-dimensional 
gridded data sets of the primary atmospheric 
variables which can be used to periodically 

United area forecast -dels While most 
national weather centers use global forecas 
model products to provide the first 
the analyses used to initialize their 
area forecast models, there is no reason why an 
equivalent data assimilation system should not 
be operated on a limited area domain provided 
care is taken in the specification of the has 

lateral boundary conditions. Hills (19S1) ft 

sv.tes. using FCCE sets ov«r . U.lt.d *”* 

ro-ain in the Australian region. 

The essential phases of data assinilat on 
scheme are a model forecast to the data 
insertion time, the correction of the forecas 
fV'-Jc to fit the observed data (the analysis 
phase) , and an initialization phase to balance 
the adjusted fields prior to the next ^eca 
crate . Lorenc (1984) lists the following 
desirable principles of such a scheme: (1) the 


analysis must fit the observations to «Uhin 
their estimated observational errors; (2) the 
analyzed fields must be internally consistent, 
matching the structure, scale, and balance of 
the atmosphere; and (3) the anlysis must be near 
the forecast based on earlier observations 
unless current observations indicate otherwise. 

Two further requirements were made. 

First, the svstem should be computations y 

md t «ogr.pl>lmU7 purokl.. s.coud, 

the svstem should be as robust as possible, 
without excessively compromising the quality o 
the output product. That is, minima man 
intervention should be needed in the operation 
of the system. This is, however, to -me exten 
at the expense of maintaining maximum detail in 
tie analysis, and a variety of options exist to 
relax these constraints in development and 
f sef Inff of this BVS t CO • 

testing^ ^ decided ^ fcase the flrgt genera- 
tion of this system on existing analysis and 
forecast modules, with which some and 

expertise was available, and to thes * . . 

pre-existing components to improves their suitabi- 
lity as assimilation vehicles. ™e basic 
modules which are used in the scheme are the 
three-dimensional variational analysis scheme 
R. S. Seaman (e.g. Mills. 1981) 

Australian Region Primitive Equation (ARPE > 

model which has been extensively descri - 

Leslie et al. (1985). The interface programs 
have been extensively modified so that the 
svstem should approach satisfying °* ea ^ 
condition (3) above. The range of «v«il»ble 
options Is summarized in Table 1. * a 

incremental assimilation cvcles are schematicai 
ly represented in Figs. 1 and 2. respectively. 
The analysis scheme, the forecast model. the 
Interface codes, and how they are linked tegeth 
er to perform an assimilation sequence in the 
McIDAS environment will be discussed in detail 

in a forthcoming paper. 

The assimilation system was ir.itlall. 
tested in a case study mode using a mixture of 
VAS soundings and special network raalosondq 
data at 1430 and 1730 GMT on March 1982 (Diak. 
et al., 1966). To adequately test the 
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Table 1 

Options To Be Selected 


OF POOR QUALiTY 


tata_jtVs : 

JTcTTscnoe 

feither/or) 

Cloud and watejr vapor drift vinds 
SVCA/buoys , 

Assimilation type : 

Basic 

Incremental 

A nalysis : , , , 

Weight ^of geostrophic/gradient wind 
Rejection tolerance 
Minimum pass radii for first pass 

Pre-processor : 

Sirma level disposition 
Divergent or non-divergenC vinds 
Telegraphy enhancement factor 

Model : 

Timescep 

Boundary conditions - fixed 


Initialization 


Physics options 


- updated 

- none 

- free VM1 

- pressure constrained VMI 

- surface T, Tp prediction 

- convective pptn 

- large scale pptn 

- vertical diffusion 



Pig. 1. Basic (or cold start) case study 
configuration. 

assimilation system using data from more 
extended periods of tine, advantage was taken or 
CTMSS support of the National Hurricane Center 
to assimilate radiosonde, buoy, VAS temperature 
and devoolnt (Smith, 1983) and gradient v n 
data (Hayden, 1985). and cloud drift and water 



Fig. 2. Incremental assimilation sequence 
configuration. 

vapor motion vinds (Stewart, et al., * 

12 hour assimilation interval run at 00 and 12 
C if r w .„ performed in a quasi-real-time mode fro 

lf<£ %S" o”ob.r 1955 

1985. No VAS or cloud and water vapor drift 
vinds were available during the weekend. For 
convenience, the domain was chosen to be th 
same as the March 6 assimilation: 26x3 grid 

with a 125 km grid spacing « n “ r * d on lon * ltude 
95W and between latitudes 20 and SON. 

A significant amount of d; ay-to-day ± 

fluctuation was noted during this per o 
«„ Inc. resting .!»»•« 

W U. « 12 OH » 4st 253. ch. 

makes verv large changes to almost all fields, 

and dearly there was c ° nslder f le f ^oied C bv the 
the initial state, which was then followed 
Ursest 1000-500mb thickness errors in any of 
12-hour forecasts (Fig. 3). It vas encouraging, 
at this stage of development of the system, to 
find that this did not destabilize the .!!!$“* in 
There is a considerable degree of oscillation in 
the skill of the 1000-5G0mb thickness £««««* 
based at 12 GMT and those based at 00 CMT. 
cause of this may be due to the “ ixes 

of satellite and conventional data at ™ 
different times, or may be due to more subtle 
factors needing further Investigation^ 



— « ... no 

a«a ‘ radiosonde data 

Fie. 3. Root-mean-square fit or racier 

1000-500 mb thickness to assimilation eye 

1985 (day 282) to 12 GMT 18 October 1985 (day 
291). 
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When looking at the RMS vector errors in 
Flc. 4, It must be remembered that the wind 
fields input to the prognosis model are non- 
divergent increments added to the model forecast 
wind fields, and thus the RMSVE statistics have 
a creater magnitude than would be expected from 
a scalar analysis of wind components. The RMSVE 
of tne radiosonde data to these wind fields are 
also included for comparison, and are substan- 
tially lower than the statistics for the winds 
used for model initiation. 



Fig. 4. Root-mean-square vector error at 500 mb 
of 12-hour prognoses, analyses with non-divergent 
wind increments, analyses with divergent wind 
increments (triangles), initialized fields, and 
12-hour prognoses at 12-hour intervals from 00 
GMT 9 October 1985 (day 292) to 12 GMT 18 October 
1985 (day 291). 

To a large extent, Lorenc’s condition (3) 
is met by correcting the first guess forecast In 
the anlysis. Of equal desirability is that 
Information from the prognosis model which is 
rot resolved by the observations or the analysis 
scheme be preserved through the analysis (correc- 
tion) and initialization phases of the assimila- 
tion. One such field is the vertical motion 
structure. The 500 mb omega field (3p/3t) were 
compared before the analysis and after initiali- 
zation at each data insertion time and, in 
general, the major features and structure were 
preserved which was very encouraging. Two 
examples for the 12-hour forecast vertical 
motion and the after initialization vertical 
motion at 00 GMT 10 October and 00 GMT 16 
October 1985 are shown in Fig. 5. In each case, 
it is clear that substantial continuity of these 
fields has been preserved, even though there has 
been significant changes to the prognosis fields 
during the analysis. These patterns are quite 
typical of the synoptic sequence during this 
period, with an area of upward vertical motion 
toward the northwest of the grid, associated 
with an advancing trough, and another weaker, 
aloncated southwest-northeast area of upward 
motion which has propagated southeastward in 
association with the eastward passage of the 
preceedine middle latitude trough. 

This period was not synopticallv very 
active. A powerful anticvclonic circulation 
covered the southeastern part of the U.S. 
throughout the entire period, with a sequence of 
‘ c*. -nclng troughs which entered the analysis 
c^rain from the vest and northwest and were 
.orced northward by the anticyclone as they 
coved eastward. The assimilation sequence 
clearly demonstrates a synoptic continuity 




Fig. 5. Before analysis and after initializa- 
tion 500 mb vertical motion (3p/3t) fields at 00 
GMT 16 October 1985 (units mb/hr). 

throughout the entire assimilation period for 
all variables, as would be expected from a robust 
analysis/prognosis system. Figure 6 shows an 
ex am ple of four assimilation analyses of 500 mb 
heights and winds for the tvo-day period 00 GMT 
15 October to 12 GMT 16 October 1985. 

Development of a functional, stable 
assimilation system is being continued in a 
delayed real-time mode in association with the 
special GALE observing networks at 12 and 6 
hour analysis/prognosis assimilation cycles, and 
in a quasi-real-time mode on a 3 hour assimila- 
tion interval during the special observing 
periods. 
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1 . 


INTRODUCTION 


Radiosonde observations near the surface 
and aloft often detect pressure, temperature, and 
wind patterns associated with large scale atmospher- 
ic wave features which govern our dally weather 
conditions. The detection of small scale waves 
such as gravity and bouyancy oscillations however, 
usually require special observations to measure the 
perturbations directly. These waves need not al- 
ways be measured or sensed directly and their ob- 
servable signatures (i.e., clouds) can be detected 
in visible photography or satellite Imagery, 

The detection of mountain and lee wave 
structures is of:en done by observing cloud sig- 
natures In regions where vertical oscillations are 
forced by terrain. Satellite observations of these 
small scale clouds are possible with high resolution 
visible sensors. Direct observation of wind and 
thermodynamic parameters is achieved by Instrument- 
ed aircraft, Ildar, and radar measurements. From 
these observations, quantitative information la 
available about mountain and lee waves, their form- 
ation and three dimensional structure. However, 
because s moist environment Is not necessary, man y 
wave features may go undetected by these convention- 
al measurements because of the absence of clouds. 

During one of the flights of the Multispec- 
tral Atmospheric Mapping Sensor (MAMS) on a NASA 
U2 high altitude aircraft in 1985, wave features 
were found in water vapor channel Imagery in close 
proximity to the Sierra Nevada mountains of Cali- 
fornia (Fig. 1). Visible channel imagery from the 
flight was examined for cloud features but none 
were found. This paper presents evidence that the 
waves in the imagery reflect changes in the vert- 
ical temperature and moisture distribution pro- 
duced by the vertical oscillations. These features 
on the downwind or lee side on the mountains are 
called lee waves. 


W. Paul Menzel 
ROAA/ASPP 
Madison, Vise. 


Robert J . Atkinson 
General Electric Co. 
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Table 1 


MAMS: The Multlspectrsl Atmospheric Mapping £ensor 

Provides high resolution visible sod 
XR scanning capabilities from a high 
altitude aircraft platform. 

Scan rate 6.25 or 12.50 rps 
Instantaneous f ie Id- of- view 2.5 or 5.0 mrad 
Crowd r solution (nadir, 19. 8km) 50m or 100m 
Total fleld-cf-vlew 66° 
Roll correction 115* 
Calibration sources IR, 2 controllable bbs 
Pixels per scan line 716 


Channel Ws 


Spectral Bands 
t length Spectral Region 


'‘tolSeflf’ 1 ' 



(microns) 


2.5/5.036.2 

1 

•42 - .45 

blue 

1.0 

2 

.45 - .52 

blue/grees 

.5 

3 

.52 - .60 

green 

.4 

4 

.60 - .62 

ye llov/o range 

.5 

5 

.63 - .69 

red 

.5 

6 

.69 - .75 

red/near IR 

.4 

7 

.76 - .90 

near IR 

.4 

6 

.90 -1.05 

near IR 

.4 

9 

6.20 - 6.90 

water vapor 

.40/. 20 

10 

10.32 - 12.02 

window (cold) 

.15/. 10 

11 

10.32 - 12.02 

window (warm) 

.15/. 10 

12 

12.20 - 12.56 

water vapor 

.75/. 40 

(1) HEAT (cu/s ter/ n*/e« ') for 1-8, 

HEAT fC) for 


2. BACKGROUND 

2.1 Instrument Characteristics 

The MAMS la a multlspectrsl scanner capable 
of producing high resolution iaagery in visible 
and infrared bands with detailed spectral resolu- 
tion. Table 1 summaries the characteristics of the 
instrument and its spectral bands. The MAMS la 
flown on a NASA U2/ER2 at an altitude of 20ka pro- 
viding a resolution cell of 50-100m on the ground. 


The width of the entire cross track field-of-vlev j 
1® 40km, providing detailed resolution over a rel- 
atively large area. j 

, Eight channels of the Instrument are used to 

separate the visible and near-inf rared radiation < 
into narrow bands which sense reflected solar energy 
from the earth's surface and clouds. The three 
Infrared channels (one of which is redundant) cover 
selected regions where infrared emission from the , 
earth and atmosphere are of interest. These bands 
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IHFIARED VISIBLE 


(■tllD-TROPOSPHERIC (CLOUD/LAND 

MOISTURE) SCENE) 



Figure 1. 


KAKS water vapor and *i*lbl* channel ccnpoalte Imagery from four adjacent flight track* on j 
22 January 1985 from 1942-2109 CHT. The lnaet ahovs the major terrain features of the 
Nevada mountains over which the imager; was collected. J 


include the mid-tropospheric (channel 9, 6.5 microns) 
and low- tropospheric (channel 12, 12.3 microns) water 
vapor region, and a window region (channel 10/11, 

11 microns) where the atmosphere is relatively trans- 
parent. The 11 and 12 micron bands fox a what has 
become known as the split window channels ( Chesters 
et al . , 1983). Further details about Che instrument 
and spectral bands say be found in Jedlovec et al. 
(1966). 

2.2 Structure of Lee Waves 

Stationary waves such as mountain or lee waves 
are a local phenomena which may have significant 
effect on local weather conditions because of their 
associated cloud patterns and wind distributions. 

They may also feedback to the larger scales of mo- 
tion. The presence of either mountain or lee waves 


I 

Is due to vertical oscillations of air parcels as' 
.air la forced over the mountain. Mountain waves 
usually consist of one vertical oscillation of the 
.flow pattern between two mountain peaks and are 
mainly governed in size by the spacing of the ridges. 
Lee waves f although forced by mountains, are not 
| restricted In size to one or two vertical oscll- • 

: latioQs as with the mountain waves. In fact, the 
overshooting of the downward flow (the restoring 
force Is gravity) on the lee side of the obstacle 
produces vertical oscillations which may extend for 
'many wavelengths downstream. 

The water vapor Image In Fig. 1 shows wave 
structures in the lee of the Sierra Nevada moun- ! 
tains. Three or four waves are dlscemable, each 
having an estimated wavelength of 15- 20 km. The 
> waves extend away from the face of the aountaln f 

♦ ran 8 c f° r about 70km and along the face for at least 

» 
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.Okjn. The vtvei appear somewhat Irregular, poa- 
lbly due to variation# In the mountain peaks, tsn- 
jental diatortion fro* the scanner geometry 
(Jedlovec et al .. 1S86), limb brightening (due to 
Increased water vapor path length at the edges of 
the flight track), as well as other unknown causes. 

The nature of the vertical oscillation pro- 
ducing lee waves was studied In detail by Brunt (1927) 
who derived a set of equations which can be u*ed to 
describe this flow pattern. Froj Brunt’s work, a 
staple vertical oscillation in the atmosphere Kna 
a period given by 



where 


(Xc-np 


% Is the envlronaental lapse rate, 

£ the dry adiabatic lapse rate, 

g the gravitational acceleration, and 

T the environmental temperature. 



An air parcel moving horizontally with a velocity 
of V end under this vertical oscillation will 
complete one cycle over a distance of ; 


)=vt = 


2iVv 


[4(r-r)] v 


( 2 ) 


Upon examining this equation. It becomes apparent 
that the horizontal wavelength (JO Is governed by 
the wind speed and atmospheric stability. Strong 
winds and low stability favor long wavelengths, 
while weak winds and strong stability favor shorter 
ones. 


K number of investigators have studied lee 
waves by observing their visible signatures and 
with direct measurements of temperature and wind 
velocity. Atkinson (1981) presents a fairly ex- 
tensive review of this research. The character- 
istics of these waves vary from occurance to oc- • 
curance but representative values are presented j 
below. The wavelengths of lee waves range from 
2-7 Okm but most occurances are In the range of 
5-20km. The wavelengths tend to increase with 
height and seem to be related to the mean wind 
speed. Corby (1957) performed s linear regression 
of wind speed on wavelength and suggested that a 
high enough correlation existed for the wavelength 
to be determined from the wind apeed alone. 

Another feature which is of interest is the 
wave amplitude, defined as the vertical distance . 
between the peak and trough of a streamline. This 
Is more difficult to determine from the cloud it- 
self and requires a number of special measurements 
of the wind field. From the available empirical 
data, these wave amplitudes are typically about 
lkn but amplitudes of 7ka have been measured. The 
larger amplitudes may occur when the wavelength 
closely matches the size of the obstacle forcing 
the wave. 

From these studies. It appears that two , 
parameters control the appearance and size of lee 
waves; namely, stability and wind. Of course the 
obstacle (mountain) is necessary to initiate the 
vertical oscillation and may have some control over 
the wave amplitude. Stability and wind are also 
important parameters in (1), however, the import- 
ance of stability in determining wavelength seems 
less important than indicated by (2). A typical 
environment al temperature and wind profile would 


consist of layers of varying stability sod an ab- 
sence of horizontal wind shear. In the lowest 
layer, the atmosphere would be neutral with a 
.ctrong layer of stability aloft. This stable lay- 
er la accessary for the formation of the waves and 
should occur near the level of the mountain peak. 

A layer of speaker stability usually exists above 
this level. The flow should be perpendicular to 
•the mountain at all levels and increase in speed 
with height. Sif-n if leant horizontal shear reduces 
the llklihcod for the development of lee waves. 

I 2.3 At mospheric Conditions j 

I In the previous section, the atmospheric 

{conditions required for the development of lee waves 
ivere suxaacrlxed . Although they are local conditions, 
j they are governed by the large scale atmospheric , 

, patterns. Thus, a discussion of the large acale - 
jfe~-rure£ on 22 January 1985 corresponding to the 
.MAKS observations is necessary. The upper air : 
:data (Fig. 2) exhibited e pronounced ridge of high 
; pressure over the eastern Pacific ocean end over 
! the northwestern portion of the U. S. A trough was 
'positioned over the central part of the country 
with e portion of it hanging back Into Utah and the 
j southern portions of Nevada and California. This 
'set up a prevailing flow pattern from the northeast 
• over much of the Sierra Nevada mountains. 



b 



Figure 2. a) 500mb ravins onde station plots and ; 
height analysis, and b) 700mb plots and streamlines 
for 23 Jan 1985 at 0000 QfT. Plots and analyses 
i are in their standard units. 
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3. DISCUSSION i 

al though Fig. 1 was briefly discussed in I 
previous sections, further discussion see** appro- 




Flgure 3. SWI-hogP^disgraM^for ^ 8 ^® d 0000 
IZ'lTtl 1985? Winds «; Plotted in knots. 
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». RADIATIVE TRANSFER SIMULATIONS 

The water vapor Image presented In Fig. 1 
represent* energy emitted by atmospheric water va- 
por somewhere in the mid-troposphere. The pro- 
cedure to determine this layer or level has been 
the subject of Investigation by researchers for 
many years. In the case of cloud height determina- 
tions using window channel radiances, the altitude 
of the cloud can be estimated by comparing the 
measured blackbody temperature to a local ravinsoode. 
The water vapor signal does not come from a single 
level but actually a layer of the atmosphere 
(Jedlovec, 1964), and depends on the amount, ver-' 
tlcpl distribution, and temperature of the vapor 
Itself. It la unclear whether blackbody temper- f 
atures In this channel can be used directly to In- 
fer a mean height from which the signal originates. 

In an attempt to estimate the amplitude of 
the lee waves observed in Fig. 1, several approaches 
were considered. Table 2 presents the observed 
blackbody temperatures measure in the peaks and ! 
troughs of the lee wave features. Typical values 
Indicate a 6* temperature difference from the peak 
to the trough of the wave. It is difficult to use 
one of the soundings in Fig. 3 to directly infer 
the levels of the peak and trough of the lee waves. 

It may be reasonable to assume that the level of 
the neutral point of the wave is 246 # K (-27*C) and 
the peak and trough oc curing at the level of their 
measured temperature (Table 2), however. If this 
is done, the estimated wave amplitude Is about 1km. 
While this Is consistent with the amplitudes of a 
typical lee wave, it depends heavily on the lapse 
rat>* of the sounding used. Because of the lifting 
of the lower level inversion by the mountains, a 
smaller lapse rate of temperature probably existed 
over the Sierras, producing a greater vertical 
amplitude for the lee wave. 

Another approach to estimate the lee wave 
amplitude used a composite of the tvo soundings in 
Fig. 3. From the composite sounding, MAMS vater 
vapor channel radiances were simulated (McMillan 
and Fleming, 1976; Fleming and McMillan, 1977; and 
Weinreb and Neuendorf fer, 1973). In order to ob- 
tain a repre entative profile in both the peak and 
trough of the waves, the Inversion layer near 700mb 
and the accompanying moisture distribution were 
raised and lowered adlabatlcally in increments of 
1 km. Blackbody temperatures were cot-puted from 


Table 2 

i bserved and simulated blackbody temperatures 
for the lee waves. Uaits are in degrees* k. 


Observed Blackbody Temperatures 

mean extreme 

Wave peak 243.5 242.0 

Wave trough 249.5 251.0 

Simulated Blackbody Temperatures (4ki;) 

Wave peak 243.9 

Wave trough 249.5 


the channel radiances. Temperature differences 
of 6*K were observed (similar to the measure values 
from Fig. 1) as well as the appropriate blackbody 
temperatures when the inversion was at 2 and 6km 
.level for the trough and peak, respectively. This 
4km vertical amplitude, while larger than most lee 
;wave amplitudes, la more consistent with the level 
of the mountain peaks which produce thea. . 

! 5. CONCLUSIONS | 

j The above work has presented evidence for j 

the existence of lee waves in the absence of clouds 
jusing high resolution MAMS moisture imagery. Al- 
though verifying measurements do not exlfct, -n eval- 
uation of wind profiles and stability in the sur- 
> rounding region support the existence of these waves. 
! Estimations of the horizontal wavelengths of the 
(observed lee wave features are hot inconsistent 1 
*with those from previous studies. The estimated . 
j vertical amplitude, el though less certain, seems 
j realistic for the setting. Further work is neces- 
sary to more eccurstely determine this vertical i 
amplitude. 
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Hurricane Debby — An Illustration of the 
Complementary nature of VAS Soundings 
and Cloud and Wa f er Vapor Motion Winds 


John F. Le Marshall' 
William L. Smith , 
and Geary M. Callan 3 


AbvJrcct 

The utility of VISSR Atmospheric Sounder (VAS) temperature and 
moisture soundings and cloud and water vapor motion winds in de- 
fining a storm and its surroundings at subsynoptic scales has been 
examined using a numerical analysis and prognosis system. It is 
shown that the VAS temperature and moisture data, which specify 
temperature and moisture well in cloud-free areas, arc comple- 
mented by cloud and water vapor motion data generated in the 
cloudy areas. The cloud and water vapor “winds” provide thermal 
gradient information for interpolating the soundings across cloudy 
regions. The toss of analysis integrity due to the reduction of VAS 
sounding density in the cloudy regions associated with synoptic ac- 
tivity is ameliorated by using cloud and water vapor motion winds. 
The improvement in numerical forecasts resulting from the addition 
of these data to the numerical analysis is recorded. 


1. Introduction 

Between 14 September 1982 and 20 September 1982, Hurri- 
cane Debby moved through the western Atlantic Ocean, 500 
miles from the east coast of the United States providing a first 
opportunity for the VAS instrument on GOES-5 to observe a 
hurricane and its environment. During this period, a data 
base was established from the usual conventional observa- 
tions, VAS observations from GOES-5, HIRS and MSU ob- 
servations from NOAA-6 and NOAA-7, cloud and water 
vapor motion “wind* 1 * observations, and National Hurricane 
Research Laboratory aircraft dropwindsonde measurements. 
The utility of the various components of this data base in de- 
fining the stor.n and its surrounding environment has been 
examined. This study concentrates on the utility of the 
GOES-5, VAS temperature and moisture sounding data, 
and cloud and water vapor motion data. The timeliness, den- 
sity, and coverage of GOES sounding and wind data make 
them well-suited for analysis and forecast purposes in the 
western Atlantic Ocean. 


2 . The data base 

The data base used in this study had three components: 
the usual conventional observations without ship reports, 


AIREPS, and the satellite components. These data were re- 
ceived in real time at the Space Science and Engineering Cen- 
ter (SSEC) at the University of Wisconsin in Madison. 

VAS sounding data from GOES-5 were processed using 
the physical retrieval scheme developed at the NESDIS De- 
velopment Laboratory in the Space Science and Engineering 
Center at Madison (Smith, 1983). This scheme provides 
soundings that represent a full physical solution of the radia- 
tive transfer equation. Surface data are used to reduce sound- 
ing errors near the lower boundary, and the temperature and 
moisture retrievals are made in a fully interactive Man- 
computer Interactive Data Access System (MclDAS) envi- 
ronment (Suomi et at., 1983), which among other things, fa- 
cilitates the removal of errors due to cloud contamination 
and noise. The first guess field used in the retrieval Scheme 
was a 12-hour National Meteorological Center (NMC) 
global forecast, and the retrievals were made approximately 
every 150 km. The quality of these VAS soundings is dis- 
cussed in detail in Velden et at. (1984) and shows that the 
RMS temperature difference between these soundings and 
radiosondes is between 1 K and 2 K through the troposphere. 

Cloud and water vapor motion winds were also produced 
at the SSEC on the MclDAS using the technique described 
by Mosher (1978), which was applied to GOES-5 half-hourly 
visible and infrared imagery and hourly water vapor im- 
agery. This technique is based on manual selection of suita- 
ble cloud and water vapor tracers from which pursuant pro- 
grams generate wind vectors, using a correlation technique 
to determine the final position of the tracers selected. The 
quality of these cloud drift and water vapor winds is dis- 
cussed in Mosher (1978) and Mosher and Stewart (1981), re- 
spectively. The height assignment for the cloud drift winds 
was done using the bi-spectral method, while the water vapor 
winds relied on the use of the water vapor brightness temper- 
ature and nearby temperature-pressure soundings. A most 
important advantage of the latter remotely sensed da. a t>pes 
was that when some areas appeared to be data-sparse during 
analysis, a quick reexamination of the VAS radiance and 
cloud and water vapor image fields often allowed generation 
of supplementary sounding or wind data where they were 
required. ~ 


1 Visiting Scientist from Bureau of Meteorology. Australia. Pres- 
ent affiliation: Cooperative Institute for Meteorological Satellite 
Studies 1225 West Dayton St., Madison. WI 53/06. 

: NOAA/NESDIS Development Laboratory, 1225 West Dayton 
St , Madison, WI 53706. 

* NOAA/NESDIS Development Laboratory, Systems Design 
and Applications Branch, 1225 West Dayton St., Madison, WI 
53706. 

£ 1985 American Meteorological Society 


3. The analysis-prognosis system 

The chief characteristics of both the analy sis and prognosis 
schemes are listed in Table 1. The analysis and prognosis 
schemes have a horizontal grid spacing of either 67 km or 1. 
km (51 X 61 or 24 X 29 elements) on a Lambert Conformal 
Conic Projection. For this study, the prognosis model mte- 
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Table i. Principle features of the analysis scheme 
. prognosis model. 


Analysis System 

- »**- 

m\'Tm ST’ » - •». w. 5 oo. «o. 300. jio, 20 o. 

Horizontal resolution: 67 km 

Ftelds analyzed or derived at each pressure level- 

• geopotential height ' 

• temperature 

• dewpoint 

• wind components 

• stream function 


Prognosis Model 

Primitive equations model in ^-coordinates 
Ten vertical levels at a = 0.05, 0 15 0 25 

Horizontal resolution (67 km) 134 k m 

daggered horizontal grid (Arakawa “C" grid) 
Semi-imphcit time differencmg 8 ' 

Br^ Ce/b I 3Undary ' ayer P arame terization 
Broad scale precipitation 

r^' t , y ^ C ° nVeCtive ^rasterization 
Updated boundary conditions 


0.95 


- 

a. The analysis scheme 

£S~=S==- 

riota ^r„T“ n "*' h ' ieh '- ,hick " c “' «"i 

E^SHSSE 

,r 1, " y " d " ,,Ckn '““ above „ir 

<*-"■ information into V’'®"' 

^r^rrr d8radi '" ,, "“ i -^ 
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The wind fields are derived from the blended geopotential 

ndtng procedure (Seaman 1977 WnH * 

fane, ion field a, eaeh analyai, [,,ei n^ J? * 

Its integration using this non-divergent field C ° mmenCeS 

b. The forecast model 

iS , a PrimitiVe eqUation model originally 
a/. 1981) The mnl. r? Amer,can rt 8 l0n (NA) (Mills et 

in a research mode within the ANMRCat 125 km Th a 
sionTnthr Cr thC NA rCg, ° n difTerS fr ° m thc Austr *Iian ver- 

1974) has CO " V ' C * Ive P ar ametenzation scheme (Kuo 1965 

schema h S i » ** (I978 ^ A comprehensive boundary layer 
cheme has also been added to the model (Mills et al 1983) 
using simtlarity theory at the surface coupled S ven ca l 


4. The a lalyses 

(C + V), and convent, onal plus VAS sounding plus cloud 
and water vapor drift wind data (C + V + CD W) The r 
guess fields for the 0000 GMT 16 Seotember loso^' f 

(NMO C 7 h b f ,,0n ° f the National Meteorologica| a Center 

LFM arfd 3nd Llm,tCd Fine Mesh (LFM) analyses The 

alobl T 7" U$ed by ,he ana,ysis sch ^‘ to modify a 

ttVedTr St ' gUfSS fie,d - The ?'obal analyses alone 
re used as first-guess fields for the other analyses due to 

problems generated by the proximity of the storm to the 
FM boundary and the better quality of the global analyst 

of h TiT! ° f ^ StOTm ‘ and because the later c ut-off rime 
other data.* inC,usion ^«*cllite and 

,J h ; m . Mn sea levcI P rcssure (MSLPi analyses for the latter 
o analyses were very slightly altered by the addition of 

cy T e — 

^rm. By using these SS fidTK^£ ^7^ 
some ways gauged the impact of VAS and CDW data used in 
MSIpTm W “ h a " “ opera ‘ ional ” MSLP analysis The 
particular rime. S "" fW SC ‘ ° f ' anaJ ^« d one at a 

Sb ?7 1* 1 ' con ventional (C), conventional 
con h sounding (C + V) and conventional plus VAS 

(C + v + P CDwrr d 5nn d T'" V3P ° r dfift Wind analyses 
) at 500 mb and 250 mb at 0000 GMT on 16 




Fig. 1. (a) Comparison of 500 mb geopotemial analyses at ^ GMT on 16 

September 1982 using conventional (C). conventional plus VAS (C + ), 

conventional plus VAS plus cloud and water vapor drift wind observations 

(C + V + CDW). (b) Comparison of 250 mb geopotential analyses at WOO . 

on 16 September 1982 using conventional (C). convent.ona plus VASjC + V) 
and conventional plus VAS plus cloud and water vapor drift wind observations 
(C + V + CDW). 


FlO. 2. (a) The 250 mb geopotential analysis at 0000 GMT 16 September 1982 
using conventional and VAS observations (C + V). Cloud drift winds and VAS 
sounding positions superimposed, (b) The 250 mb geopoteni.al analysis at 0TO0 
GMT on 16 September 1982 using conventional. VAS ..and cloud drift W'^obser. 
vations (C + V + CDW). Cloud drift winds and VAS sounding positions 

superimposed. 
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Fig. 3. (a)The250mbgeopotentialanalysisat0000GMTI7September 1982 

using conventional and VAS observations (C + V). Cloud drift winds and VAS 
sounding positions superimposed, (b) The 250 mb geopotential analysis at 0000 
GMT 1 7 September 1982 using conventional, VAS, and cloud drift wind obser- 
vations (C + V + CDW). Cloud drift winds and VAS sounding positions 
supenmposed. 
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September 1982. They illustrate the impact of increasing the 
data base and show a gradual correction of the geopotential 
fields to what on - would expect to be an optimum analysis 
with all data being used. Figs 2(a) and (b) show the 250 mb 
C 4- VandC + V-f CDW analyses for 0000 GMT in a little 
more detail, with the VAS sounding positions (+) and cloud 
and water vapor motion wind data plotted. Fig. 2 illustrates 
two salient features of these analyses. The first is that the 
VAS sounding distribution and the cloud and water vapor 
motion wind distribution are complementary and provide a 
quite uniform data base over the western Atlantic. This is an- 
ticipated, as the VAS sounder provides soundings in the 
cloud-free areas, while suitable cloud tracers can often be 
found in the remaining areas. The second is that although the 
VAS data improve on the conventional analysis, the loss of 
VAS data caused by the intense cloudiness associated with 
Hurricane Debby has masked two significant features in the 
analysis (Fig. 2(a)): the exact positions of the cyclonic circu- 
lation centers near 70°W, and the strength of the ridge at 
60° W. It can be seen that the addition of cloud drift wind data 
to this analysis (Fig. 2(b)) has ameliorated the problem. Sim- 
ilarly, in the 0000 GMT analyses on 1 7 September 1982 at 250 
mb, the VAS sounding and cloud drift wind distribution was 
again uniform (Figs. 3(a) and (b)), while the addition of 
cloud drift wind data again ameliorated the problem of defin- 
ing ridge and trough positions in an otherwise data-sparse 
region. It should also be pointed out that during the 0000 
GMT analysis on 17 September 1982, additional cloud drift 
wind data were generated from imagery in a matter of min- 
utes after it was found a data void was present near 38°N 
67° YV. 


Table 3. SI skill scores of the 12-hour forecasts delineated 
according to the data base of the original analyses. 


Level C C 4 V 


C 4 V 4 CDW 


MSLP 48.6 46.4 

1000 49.9 47.6 

500 49.9 49 j 

250 51.1 43.8 


46.2 

46.5 

49.8 

41.7 


tude of these position errors, but the relative magnitude gives 
an indication of the positive impact of both VAS sounding 
and cloud and water vapor wind data on the analyses. The 
MSLP 12-hour forecast fields from the CandC + V4- CDW 
analyses at 1 200 GMT 16 September 1982 are shown in Figs 
4(a) and (b). 

SI skill scores (Tewcles and Wobus, 1954) where SI = 
100£e f /Zg«, e g is the error in the forecast variable gradient, 
and g e is the observed or forecast variable gradient (which- 
ever is larger) have also been calculated for each of the model 
integrations and are presented in Table 3. When examining 
these scores, it appears the difference in skill shown in the po- 
sitioning of the center of the storm is partially reflected in the 
SI scores at MSLP and 1000 mb. It must be noted, however, 
that the surface analyses have been constructed without sur- 
face reports over much of the analysis area, hence the com- 
parison of gradients in these regions must be viewed with 
some degree of caution. It is significant that the 250 mb level 
is the level that contains the most data, and forecasts at this 
level have shown a steady increase in skill as the data base 
increases. 


5. The forecasts 

The analyses using C, C 4- V, and C4V4 CDW data at 
0000 GMT and 1200 GMT on 16 September 1982 have been 
used to initialize the primitive equation model, described in 
Section 3b, over the analysis domain. The forecasts have 
been run at 137 km resolution. The integrations were run for 
a period of 12 hours, were non-nested, and the C + V 4* CDW 
analyses were used for verification. It should be noted that 
with analyses and forecasts at this resolution, we can exam- 
ine only the larger-scale features associated with the hurricane. 

The central position of the storm system has been meas- 
ured from the 12-hour MSLP forecast and compared to the 
estimated best track position. The results are shown in Table 
2. No great significance is intended to be fixed to the magni- 

Table 2. Storm central position error compared to 
estimated best track position for the 12-hour forecasts. 


ERRORS (N.M.) 

Forecast valid C C4V C + V + CDW 


I200GMT 

16 Sept 1982 

91 

30 

28 

0000 GMT 

17 Sept 1983 

59 

23 

11 

Mean Error (N.M.) 

75 

27 

20 


6. Concluding remarks 

From this example, it can be seen that VAS and CDW obser- 
vations form complementary data bases. In areas of exten- 
sive cloudiness, CDW data allowed full subsynoptic scale 
analysis at cloud-tracer level by adding gradient information 
in areas devoid of VAS data. This often significantly changed 
the numerical analysis that had previously been interpolated 
through the area. It can also be seen that the use of VAS 
sounding and cloud and w'atcr vapor wind data, in conjunc- 
tion with an “opera tional" MSLP analysis, appeared to have 
a positive impact on 12-hour numerical forecasts, although 
it is recognized that two sets of forecasts do not provide a sta- 
tistically significant conclusion. An important point n c ved in 
this study is that through the use of McIDAS, wind vectors 
and VAS retrievals could be generated in association with the 
analysis process to allow selective augmenting of the analysis 
data base in a time frame of minutes when data voids were 
encountered. This capability is one most operational ana- 
lysts would find very useful. 

As a result of the analysis improvements seen in this study, 
full resolution model integrations (67 km) are now being un- 
dertaken on a number of tropical storms in order to further 
quantify the impact of VAS sounding and cloud and water 
vapor motion data in a variety of cases. Similar studies have 
already been done over the mainland, where careful account 
needs to be made of the relative quality of VAS temperature 
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gradients and absolute values. Again the complementary na- 
ture of the VAS soundings and cloud and water vapor drift 
“wind*’ data has been evident, with the cloud and water 
vapor motion wind data providing the information needed 
for the positioning of trough and ridge lines and circulation 
centers in cloud areas devoid of VAS soundings. This has 
permitted detailed subsynoptic scale analysis between the 
more widely spaced soundings of the conventional network. 

Finally, it should be noted that the study was done with a 
data base available in “real time" when the VAS instrument 
is operational. Thus, it would appear that the prospects of 
a quasi-continuous 100-km resolution upper air data base at 
one- to three-hour frequency with a timeliness suitable for 
operational weather forecasting in the North American re- 
gion are excellent. 
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1 . INTRODUCTION 

The Multispectral Atmospheric Mapping Sensor 
(MAMS) is a new instrument designed to produce high 
resolution imagery in eight visible and three 
infrared spectral bands from an aircraft platform* 
Thermal emission from the earth’s surface, clouds, 
and atmospheric water vapor is measured at 12.3, 

11.2, and 6.5 microns at up to 50 meter horizontal 
resolution. Similar infrared spectral bands are part 
of the VISSR Atmospheric Sounder (VAS) which is 
capable of 6.9 km horizontal resolution from its 
geostationary orbit. Sequences of 6.7 micron VAS 
images have been used to track water vapor structures 
and co infer atmospheric motion vectors (Stewart et 
al., 1985). They have also been used to investigate 
small scale atmospheric moisture variability 
(Jedlovec, 1984). The high resolution MAMS data give 
credence to the small scale variations in the VAS 
water vapor imagery. 

The purpose of this paper is two-fold. We 
shall (1) describe the MAMS instrument and the 
radiance characteristics, and (2) intercompare some 
VAS and MAMS radiances and images. 


2. MAMS INSTRUMENT CHARACTERISTICS 

The MAMS is a modification of NASA’s Airborne 
Thematic Mapper (Daedalus, 1982). The MAMS has two 
major components, the scan head and the spectrometer 
(Fig. 1). The scan head consists of the primary 
collecting telescope, a rotating scan mirror, a motor 
encoder assembly, and two controlled thermal sources. 
The field of view is defined by the aperture which is 
available in two different sizes, 2.5 mrad and 5.0 
mrad. The spectrometer consists of optical elements 
which spectrally separate the polychromatic input 
energy , lenses which focus the separated energy onto 
sensors, detectors which convert optical energy to an 
electrical signal, and pre— amplif iers which condition 
the signal. The energy collected by the spectrometer 
is separated into four distinct optical paths using 
dichroic filters; three paths for the infrared data 
(one for each spectral band) and one for the visible 


and near infrared data. The spectrometer 
detectors are silicon for the visible radiation 
and mercury cadmium telluride for the thermal 
infrared radiation. The visible detector is an 
eight element array which converts optical 
energy to electrical energy. The spectral 
response of each band is determined by the 
dispersion characteristics of the prism placed 
in the optical path and the location of the 
arrav in the dispersed beam. The thermal 
infrared detectors are each housed in a vacuum 
dewar which contains a cooled long wave filter 
which defines the particular spectral band. 

Each infrared detector has an associated pre- 
amplifier to process the electrical signal. The 
spectral bands are summarized in Table 1. 

The spectral response functions for the three 
infrared charnels of the MAMS (Fig. 2) were 
specified to be similar to the 6.7, 11.2, and 
12.7 micron channels on VAS (also shown in Fig. 
2 ). 

The MAMS produces high resolution imagery 
in these visible and infrared spectral bands 
from a high altitude aircraft at 20 km as shown 
in Fig. 3. The horizontal ground resolution, 

GR, of each individual field of view is 100 

Table 1 


Rangle of MAMS Spectra) Bands 


8 vis 

ihle 

3 infrared 




(microns) 

(microns) 




.42- 

.45 

6.20- 6.85 

(uppe r 

level 

water vapor) 

.45- 

.52 

10.32-12.02 

(window") 


.52- 

.60 

12.20-12.56 

(lower 

level 

water vapor) 

.60- 

.62 





.63- 

.69 





.69- 

.75 





.76- 

.90 





.90-1 

.05 
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THERMAL THERMAL 



Fig. 1. Schematic of MAKS spectrometer and scan head. 



P 



VAS 



Fig. 2. Spectral response curves indicating the 
sensitivity of the MAMS and VAS instruments in 
the 12 micron, 11 micron, and 6 micron regions. 


meters (50 meters if the smaller 2.5 mrad optics 
are used). The total field of view for one scan 
is 86° which covers the width, w, of the entire 
cross path of roughly 40 km. The scan rate is 
6.25 revolutions per second (12.50 when the 2.5 
mrad optics are used) and produces 671 overlap, 
1-P, of each scan line with the previous one for 
nominal aircraft speed (740 km/hr). On each 
scan line 716 samples are taken so that there is 
522 overlap of neighboring pixels along a scan 
line (4% when the 2.5 mrad optics are used). 
These overlaps can be used to good advantage to 
reduce the noise in the scene by spatial 
averaging. 

Radiometric calibration of the infrared 
data is accomplished bv viewing two blackbodv 
sources and by assuming the detector response is 
linear with respect to radiance (a fractional 
nonlinearity of less than 5x10 L is observed). 
Calibration occurs every scan line by viewing a 
vara and then a cold blackbodv cf known 
temperature. For a given spectral band, the 
radiance is determined from the temperature 
through the convolution of the normalized 
spectral response and the Planck function. The 
calibrated radiances are transformed back into 
temperatures using the inverse Flanck relation 
with the nominal spectral hand center vave- 
1 ength. 

Prior to flight, the gain and offset for 
each spectral band are adjusted electronically 
in order to minimize the likelihood of satura- 
tion over hot scenes and tc maximize the effec- 
tive operating range. The 6.5 micron band is 
adjusted to cover 210 to 280°K, while the 11.2 
and 12.3 micron bands operate from 220 to 340 3 K. 
An eight hit digitizer is used so that for the 
6.5 micror band the gain is approximately .3 C K 
per count, and for the remaining surface viewing 
bands it is approximately . 5°K per count. 
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Fig. 3a. A simplified schematic of the inflight 
configuration of the optical elements of the 


MAMS. 



Truncation errors arc assumed to be roughly half 
these values. 

The infrared detectors on MAMS exhibit 
some low frequency noise so that samples taken 
at 52 microsecond intervals for the spin rate of 
6.25 revolutions per second (26 microsecond 
sampling occurs at 12.5 rps) experience noise 
contributions with significant correlation. An 
analvsis of covariance (Menzel, 1980) of the 
noise with the 5 mrad aperture reveals that 
after averaging ten consecutive samples noise is 
reduced by a factor of .54, .52, and .37, 
respectively, in the three infrared spectral 
bands in order of increasing wavelength i— 

(uncorrelated noise would have a factor of I/^IU 
or 3 ) 

Inflight single sample noise values were 
determined from multiple samples over a uniform 
portion of the ocean. For a sample size of 1UU, 
the single sample noise values with the 5 mrad 


aperture were found to be .2, . and .4°K, 
resb ective ^ » f° r the s P e ctral bands in order of 
increasing wavelength (for the 2.5 mrad aperture 
the values were .5* .2, and ,9 K) . 


3. DATA INTEF.COMPARISONS 

Fadiances for the MAMS and VAS spectral 
responses were simulated for each band and 
sensor in a radiative transfer calculation. A 
representative atmospheric transmittance was 
computed using a band model and radiosonde 
temperature and moisture profiles. Table 2 
shows that the MAMS detects radiances within one 
to two mW/ster/nr/cm 1 of the VAS detected 
radiances. The MAMS spectral responses are 
shifted to shorter wavelengths (as shown in Fig. 
2), hence less radiance is detected. This 
translates into brightness temperature 
differences of less than a degree Kelvin. 

Table 2 


Inter comparison of Simulated Data 



6 micron 

11 micron 

12 micron 

R(MAMS) 

7.86 

96.5 


109.5 

T(KAMS) 

254.8 

287.3 


286.0 

R(VAS) 

8.73 

97.5 


111.6 

T(VAS) 

253.9 

287.3 


285.4 

AT 0.9°K 

MAMS- VAS 

0.0 


0 . 6 

radiances 

in mW/ster 

/m 2 /cm 

temperatures in 


Radiances from the MAMS and VAS were 
collected simultaneously over a relatively 
isothermal area of ocean off the California 
coast on 18 May 1985. Table 3 presents the 
comparison. After compensating for the inherent 
spectral differences of the two instruments, the 
MAMS data is within .5, 2.5, and 1.5°K of the 
VAS data, respectively, for the spectral bands 
in order of increasing wavelength. The observed 
split window channel differences (11 micron less 
12 micron) are comparable, 2.5°K for MAMS and 
2.1°K for VAS and thus showing similar low level 
moisture concentrations. The observed water 
vapor channel data (6 micron) show MAMS 
detecting upper level moisture concentrations 
very much like VAS, with MAMS showing somewhat 
less attenuation than VAS. These comparisons 

Table 3 

Inter comparison of Observed Data 


6 micron 

1 1 micron 

12 micron 

R(MAMS) 6.01 

99.9 


111.1 

T (MAMS) 247.1 

289.5 


287.0 

111.1 

R(V’AS) 6.80 

97.0 


T(VAS) 246.5 

287.1 


285.0 

at mams-vas ‘ 6 K 

2.4 


2.0 

radiances in mW/ster/m 

2 / cm 1 , 

temperatures in °K 
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are quite reasonable, since the spatial resolu- 
tion of each sensor is very different (VAS at 7 
Vm and HAMS at .05 km at nadir), the navigation 
of the MAMS data is from visual siting of 
landmarks only, and non-nadir viewing of the 
MAMS has not yet been taken into account. 

Imagery from the MAMS and VAS were 
compared from flights over Oklahoma and Kansas 
on 12 May 1985. The goal of the intercoraparison 
is to verify small scale features observed in 
the VAS water vapor images with the higher 
resolution MAMS data. These water vapor inhomo- 
gereities in the VAS data have been used to 
obtain motion vectors when several images are 
considered in sequence. While the imagery is 
often too amorphous to track winds reliably with 
a correlation technique, it has produced good 
wind fields by single pixel tracking (Stewart et 
al., 1985). The MAMS data tends to verify that 
snail scale gradients (down to several single 
pixels) in the VAS data are real. Fig. 4a shows 
the VAS 6.7 micron images with the MAMS flight 
track superimposed while Fig. 4b shows the 6.3 
micron MAMS data (at 400m resolution for dis- 
play) gathered in flight from G to H. These 
images were enhanced to delineate the moisture 
features. Brighter (colder) regions represent 
clouds or relatively high water vapor content 
and darker (warmer) areas represent drier air in 
the middle layer of the atmosphere. The MAMS 
images are not limb corrected and, therefore, 
images appear brighter (colder) toward the 
edges. The moisture features compare very 
nicely; starting at G, MAMS shows a dry tongue 
giving way to a moist ridge which becomes less 
moist at H. Fig. 5a shows single pixel resolu- 
tion of the VAS image covering the MAMS track 
from E to F of Fig. 4a; Fig. 5b displays the 
MAMS data covering the boxed region of Fig. 5a. 
The temperatures from brightest to darkest 
pixels within the boxed-in area differ by 5°K 
indicating that these fluctuations are not 
instrument noise, but real atmospheric features. 
Six of the VAS features are corroborated nicely 
by the MAMS image. This lends credence to the 
assertion that the VAS distinguishes moisture 
features which are small enough to represent a 
wind vector, yet large enough to exist in 
several successive images. The features we are 
examining were discernible for at least two 
hours in the VAS images. 


The MAMS has been developed to combine 
high resolution visible and infrared measure- 
ments to study small scale atmospheric moisture 
variability, to monitor and classify clouds, and 
to investigate the role of surface characteris- 
tics in the production of clouds, precipitation, 
and severe storms. As part of the Cooperative 
Huntsville Meteorological Experiment (COHMEX), 
the MAMS will be flying this summer to gather 
data for these studies. 
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4. CONCLUSION 

The MAMS measures thermal emission from 
the earth atmosphere system in eight visible and 
three infrared spectral bands (12.3, 11.2, and 
6.5 microns) at up to 50m horizontal resolution. 
These infrared bands are very similar to three 
of the VAS infrared spectral bands. The instru- 
ment was flown aboard a NASA U-2 in Mav-June 
1985 and corresponding VAS data was collected. 
MAMS radiometric performance is comparable to 
that of VAS. The MAMS exhibits somewhat less 
attenuation from water vapor than VAS because 
its spectral bands are shifted to shorter wave- 
lengths away from the absorption band center. 
Intercomparisons of these images give credence 
to small scale variations In the VAS water vapor 
imagery. Attempts are underway to identify the 
scales of which these features must be tracked 
to depict the preconvec t ive environment. 
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Aa. Enhanced VAS wafer vapor Image with the MAMS flight '«<* 11!” Iw'cMT. 

Ab. Corresponding MAMS 6.3 micron enhanced Image gathered In flight from G to H on 18 May 
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1 . INTRODUCTION 

Over the years, a number of regional 
numerical weather prediction models have been 
developed in an attempt to forecast not only 
svnopcic scale events, but the mesoscale 
phenomena that are associated with or even 
independent of synoptic scale svstens. One 
problem with regional models, however. Involves 
the forecast of precipitation, which is a 
mesoscale event itself. The onset of precipi- 
tation in model forecasts commonly lags behind 
the observations, even when precipitation is 
observed at the initial time. To understand 
the nature of precipitation forecast deficien- 
cies, two key related variables, namely mois- 
ture and vertical motion, their initialization 
ir. a model, and their evolution during the 
forecast must be examined. In this study, 
model simulations with modifications of the 
initial specifications of moisture and vertical 
motion will be examined to determine the 
relative impact of the mesoscale component of 
each parameter on the subsequent quantative 
precipitation forecast. 

2. MESOSCALE MODEL 

2. 1 Analysis Svstea 

The analysis system and the forecast 
model have a grid of 51 east— vest rows and 61 
north-south columns on a Lambert conformal 
projection, with a horizontal grid increment of 
67 km at the standard latitudes of 20 and 50 
decrees north latitude. The domain is shown in 
Fig. 2. The principal characteristics of the 
analysis system are found in Table 1. 

The objective analysis scheme is a 
version of the analysis system which is used by 
the Australian Numerical Meteorology Research 
Centre (A.VKRC). This research analysis sv stem 
w<jS developed from the analysis scheme used 
operationally by Australia's Bureau of Heteo- 
ro.ocy (Seaman et al., 1977). The scheme is a 
combination of the successive correction method 
(SCM) of Cressoan (1959) and the variational 


Table 1 

Principal Characteristics of the Analysis 
System and Model 

Analysis System 

Combination of successive correction method and 
variational blending in three dimensions 

Ten pressure levels, p-1000, 850, 700, 500, 

400, 300, 250, 200, 150, and 100 mb 

Horizontal resolution: 67 km 

Fields analyzed w derived at each pressure 
level: 

- geopotential height 

- temperature 

- dewpoint 

- horizontal wind components 

- streamf unction 

Forecast Model 

Primitive equations model in o-coordinates 

Ten vertical levels at c-0.09, 0.19, 0.29 

0.99 

Horizontal resolution: 67 km 

Staggered horizontal grid (Arakawa ,f C" grid) 
Semi-implicit time differencing (At-5 min) 
Similarity theory surface layer 
Stability-dependent vertical diffusion of 
momentum, heat, moisture above surface layer 
through depth of boundary layer 
Surface shortwave and longwave flux modified by 
cloudiness 

Surface energy balance equation 
Large-scale precipitation 
Kuo-type convective parameteriza t ion 
Updated boundary conditions 


blending technicue of Sasaki (1958, 1970). The 
research- version contains three-dimensional 
variational blending of geopotent ial . Mills 
(198!) presents a qualitative discussion of the 
analysis sequence. 
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Fig. 1. Schematic summary of experimental 
design. The duration for each simulation is 
indicated by horizontal arrows. 

2 . 2 Forecast Model 


The forecast model was developed from 
the Australian Region Primitive Equations 
(ARPE) model (McGregor et al. v 1978). The — 
subsyr.optic Scale Model (SSH) of the Space 
Science and Engineering Center at the Univer- 
sity of Wisccnsin-Madison used in this study is 
an adaptation of the ARPE model to the North 
American region. In this adaptation on Kuo- 
tvpe convective parameterization (Kuo, 1965, 
197—) and a comprehensive planetary boundary 
laver scheme were incorporated Into the model. 
Also, a finite differencing scheme discussed by 
Cerbv et al. (1972) which lessens truncation 
effects in regions of steep topography has been 
added . Further characteristics are found in 
Table 1. 

3. EXPERIMENTAL DESIGN 

Four simulations are made with the SSM 
to determine the impact of the initial moisture 
ard vertical motion fields on the evolution of 
the precipitation forecasts. These are a 
control forecast 12-hour simulation that begins 
at 1200 CKT 6 March 1982 and three experiment 
simulations with modifications to the moisture 
ar.d vertical motion fields Incorporated at 1800 
GMT. The experiment forecasts from 1800 GMT 
are compared to the second half of the control 
forecast. The experimental design is given in 
Fig. 1. 


The 12-hour control simulation begins at 
1200 CMT 6 March 1982 and develops precipita- 
tion by 1300 GMT. The model is still producing 
precipi tat ion at 1800 GMT the time a data set 
is extracted for use in the three experiments. 



a) unsmoothed, b) smoothed. Contour interval is 
i g kg -1 . 

Beginning at this time, changes in the precipi- 
tation rates and accumulations due to changes 
in the initial specifications will be most 
accentuated. The six-hour control forecast 
from 1600 GMT 6 March, to 0000 GMT 7 March 1972 
provides the "ideal forecast" — one that contains 
mesoscale moisture and vertical motion informa- 
tion from beginning to end. 

Experiment 1 is run like the control 
forecast, except that at 1800 GMT, smoothed 
moisture fields replace the mesoscale moisture 
fields in the mode). The moisture fields are 
smoothed to make them representative of fields 
that would be obtained from analyzed radiosonde 
data, usually considered to be at the synoptic 
scale. Fig. 2 shows the original and the 
smoothed mixed ratio fields for the sigma level 
of 0.69. Experiment 1 provides the precipita- 
tion forecast for the initial state with onlv a 
synoptic scale moisture field. 

Experiment 2 is like the control forecast 
except that the initial wind fields are con- 
strained to be nondivergent at 1800 GMT. The 
nondivergent wind components are derived from 
the model winds extracted from the control fore- 
cast data set and are inserted into the model. 
All of the other variables remain unchanged. 
Because the winds are nondivergent, the initial 
vertical motion is suppressed, particularly away 
from areas of significant topography (Fig. 3). 
Experiment 2 shows the effect of a nondivergent 
initialization of winds on the subsequent 
quantitative precipitation forecast. 
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Fig. 3. Vertical notion d at 1800 GMT for 
sigma level of 0.49; a) control, b) experiment 
2. Negative values (upward notion) are solid. 
Contour interval is 2X10 6 s 

Experiment 3 is a combination of Experi- 
ments 1 and 2, containing both the initial 
smooched moisture and nondivergeot vind fields. 
This experiment is the most representative of a 
forecast beginning with synoptic analyses. 

It should be noted, however* that mesoscale and 
synoptic scale information from the six-hour 
control forecast remained in the height and 
temperature fields. 

4. CONTE CL CASE OF 6 MARCH 1982 

The purpose of this research to study the 
change? made to a precipitation forecast due to 
modifications of the initial moisture and 
vertical motion fields is achieved by comparing 
simulation of the three experiments to the 
control forecast. The performance of the 
control itself is discussed very briefly here, 
although its accuracy in terms of the observa- 
tions is not central to the analysis. 

The model performed well at 500-ob level, 
as shown in Fig. 4. The forecasted trough over 
the central U.S. is slightly behind that of the 
aralvsis, but the intensity is simulated veil. 
The forecast, however, does not capture the 
short wave trough observed west of Florida. The 
forecast of mean sea level pressure was 
excellent, with accurate placement and strength 
oi the features (not shown). 

The precipitation forecast from the model 
was not as satisfactory as the height and mean 
sea level pressure forecasts. The model's 
placement of copious rainfall east of Virginia 



Fig. 4. 500-mb geopotential heights (m) and 

500-trt relative vorticity (contour interval 
40X10 6 s negative values dashed) at 0000 
GMT; a) observed, b) control. 

and North Carolina seems accurate, based upon 
comparing Figs. 5 and 6 with radar maps (not 
shown). Precipitation output from the model is 
shown in hundredths of centimeters. The 
development of the precipitation area over Ohio 
and Pennsylvania lagged four hours behind the 
observations, yet the model greatly over- 
forecasted the amounts in this region, producing 
as ouch as 2.00 cm (0.75 in.) while the 
observations were all less than 0.25 in. The 
model's performance in the Gulf region compared 
much less favorably to the observations than did 
the forecast off the East Coast and over the 
northeastern United States. The precipitation 
band observed over Louisiana, Mississippi, and 
Alabama was not forecasted at all. No rainfall 
greater than 0.25 in. was forecasted over the 
land. Instead, the model generated precipita- 
tion out ever the Gulf of Mexico, but not until 
1700 GMT. The validity of this feature is 
questionable because of the lack of verifica- 
tion data; however, satellite imagery did 
suggest that precipitation had developed here 
earlier before moving eastward. The forecasted 
feature, however, remained stationary. 

It is unfortunate that the precipitation 
forecast produced by the control simulation was 
not more accurate. The purpose of this study, 
however, is to discuss changes in the forecasts 
brought about by rodif ications to the initial 
conditions. It is expected that the compari- 
sons between the control and the experiments 
will be meaningful, regardless of the simula- 
tion deficiencies in the control forecast. 
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EXPERIMENTAL RESULTS ANT COMPARISONS 





Fig. 6, Six-hour precipitation accumulations 
(hundredths of eta) from the control forecast 
ending at (a) 1600, (b) 0000 GMT. 


5 . 


The six-hour precipitation forecasts 
from 1800 GMT to 0000 GMT of the control and 
the experiments are compared. To facilitate 
this comoarison, three rectangular domain 
’'boxes" encompassing the major precipitation 
areas of the control forecast were chosen for 
analysis (Fig. 7). Area-averaged accumulations, 
plotted in Fig. 8, were calculated in each box 
by Burning the total accumulations and dividing 



Fig. 7. Areas of precipitation for statistical 
calculations. 


In Experiment 1* accumulations were ouch 
lower than those of the control, especially in 
boxes 2 and 3, where moisture values at 1800 
QfT had been reduced by the smoothing process. 

In these two boxes, accumulations were only 
24-30Z of the control values. One precipitation 
maximum was reduced by 501, and two other maxima 
found in the control were not forecasted in 
Experiment 1. Precipitation rates lagged behind 
the control forecast, and areal coverage was 
reduced to 68—772 of the control. 

The reduced vertical motions of 
Experiment 2 resulted in lower accumulations 
that were 40-452 of the control values for the 
three boxes. Precipitation maxima for this 
experiment were positioned veil when compared to 
the control, but most maxima were much lower. 
Area coverage was also reduced, while 
precipitation rates generally lagged behind the 
control . 

As would be expected. Experiment 3, both 
lacking initial vertical motion and having 
lower initial moisture values than the control 
forecast in the precipitation regions, produced 
the lowest precipitation accumulations of the 
three experiments. Fig. 8 reveals that Experi- 
ment 3 yielded much lower average accumulations 
and precipitation rates than the control, with 
Box 2 having the largest differences. Average 
accumulations at the end of the forecast were 
onlv 4-222 of those from the control forecast, 
and precipitation covered 6-432 of the control's 
area of precipitation. This experiment missed 
completely the precipitation centers over the 
Atlantic and in the Gulf that were forecasted by 
the control (Fig. 6b). After beginning with 
primarily "synoptic scale" conditions, the model 
was unable to recover in the six-hour period and 
regenerate the mesoscale fields necessary to 
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-eid detailed precipitation features and higher 
levels of precipitation accumulation* 

An interexperimental precipitation 
c^roarison supports the results previously 
discussed. Comparing the precipitation fore- 
casts between Experiments 2 and 3 shows the 
effects of smoothing the Initial moisture 
fields* .lust as done by the comparison between 
Experiment 1 and the control. Fig. R shows 
that the average accumulations in Experiment 3 
were, at most, 50Z of the accumulations in 
Experiment 2. Areas of precipitation were also 
reduced in Experiment 3. Comparing Experi- 
ment 3 with Experiment 1 also showed the 
effects of suppressing the initial vertical 
velocities. Again, precipitation 
accumulations, rates, and areal coverage were 
reduced significantly. These comparisons are 
indeed consistent with the comparisons between 
the control and each experiment. 

6. CONCLUSIONS 

Four simulations were made with a 
numerical mesoscale model to determine the 
effects of suppressing the initial mesoscale 
information in the moisture and wind fields on 
the precipitation forecasts. The latter half of 
a 12-hour control forecast represented an ideal 
forecast that was "ini tialiaed" with fields 
containing mesoscale as well as synoptic scale 
structure. The three experiments were degraded 
versions of this ideal initialization and their 
forecasts were compared to the control forecast 
during the period from 1800 GMT 6 March to 0^00 
GMT 7 March 1982. 

The characteristics of the precipitation 
forecasts of the experiments are not different 
enough to conclude that either mesoscale 
moisture or mesoscale vertical velocity at the 
initial time is more important for producing a 
forecast closer to that of the control. For 
example, at the end of the forecast, the 
average precipitation accumulation over a land 
area in the smoothed moisture experiment was 
67Z of the control forecast, while the non- 
divergent experiment produced accumulations that 
were 45Z of the control. Over the oceans, 
however, the smooched moisture experiment 
produced 24-30Z of the control forecast, less 
than the nondivergent experiment, which produced 
4C-44Z of Che control forecast. No general 
statements can be made about this study, there- 
fore, except to sav that degrading one or both 
cf these initial fields produces a drier fore- 
cast when compared to the control forecast. 

First cf all, comparisons were made for only a 
si>.-hour period, and this may not have been 
sufficient tine for the model to recover from 
the modifications made. Secondly, these experi- 
ments were performed on one case only, with Just 
one model, and as such, generalizations regard- 
ing the relative worth of initial moisture and 
vertical velocities are inappropriate. 

This research points out the need for 
proper initialization of mesoscale moisture and 
vertical motion information so that more accu- 
rate precipitation forecasts nay be obtained. 
Improvement of the data base that is analyzed 
should be very helpful. More detailed moisture 



for (a) Bor 1, (b) Box 2, (c) Box 3. 
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information la available from high-resolution 
satellite data, satellite cloud and water vapor 
drift winds would improve the horizontal winds, 
and thus the vertical motion derived from the 
convergence/divergence of these winds. Consider- 
able research will be required tc determine how 
to use a satellite mesoscale data to make 
short-term precipitation forecaat6 more accurate 
and timely. 
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High-Order Low-Pass Implicit Tangent Filters for use in Finite Area Calculations 


William H. Raymond 
Abstract 


High-order implicit tangent filters are developed. The implicit tangent filters possess 
highly selective amplitude response function and they can be applied relatively close to a 
boundary. Comparisons are made between the implicit and the traditional or explicit 
filters. Numerical simulations are performed to test the response in a limited area model. 
An algorithm to maximize computing efficiency is presented. All tests indicate the 
desirability and utility of the implicit tangent filter. 



